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Quelaues problèmes de mécanique quantique 


MESDAMES, MESSIEURS, 


Je regarde comme un grand honneur l'invitation de prendre la pa- 
role à Paris, la ville classique de la recherche physique, et j’en remercie 
de tout cœur mes collègues de Paris. Je ne l’ai acceptée qu'avec hésita- 
tion, parce que je ne parle qu'imparfaitement la langue française, et 
je vous demande en cela toute votre indulgence. 

Le choix du sujet m'a causé quelques difficultés. En effet, j'ai publié, 
avec mon ami et collègue Jordan, un livre «La Mécanique élémentaire 
des Quanta » dont la composition nous a pris plus de deux années. 
Je n’ai pas eu le temps d'aborder un nouveau travail car ce livre n’a 
été fini qu'au début de décembre dernier. Je me permettrai donc de 
parler de quelques chapitres de ce livre, de mes propres travaux, qui 
ne sont pas encore terminés et des recherches de mes élèves. 

Aujourd’hui je voudrais dire quelques mots de la conception statis- 
tique de la mécanique des quanta. 

Vous savez tous que la merveilleuse découverte de M. Louis de BRoO- 
GLIE sur la nature ondulatoire de la matière a donné lieu à une révo- 
lution dans le monde 'de la mécanique. La mécanique ondulatoire 
de SCHROEDINGER se rattache d’une façon directe aux idées de 
M. de BROGLIE. Tout d’abord, les deux savants ont beaucoup espéré 
qu'il serait possible d'éliminer tout à fait la conception du 
mouvement des particules (atomes, électrons, protons ou quanta 


HE 
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de la lumière) et de lui substituer des phénomènes ondulatoires. Mais 
il est démontré que cela ne peut se passer ainsi. Chaque phénomène 
peut s’interpréter de deux façons, l’une d'elles s’expliquant par la 
mécanique des particules, l’autre par la mécanique des ondes. 

Un exemple nous permettra de saisir la différence entre ces deux 
interprétations ; considérons la répartition des particules dans l’espace 
après choc contre un obstacle quelconque. Si par exemple un faisceau 
d'électrons traverse une fente ou rencontre des atomes, il se produit 
un phénomène de diffraction dans la théorie des ondes, et un chan- 
gement d’orbite dans la théorie des particules ; tous deux corres- 
pondent à des faits observables: d’une part on a pu constater des 
phénomènes d’interférence, et d'autre part on a pu aussi montrer la 
nature atomique des particules par leur dénombrement (scintillations 
ou compteur de GEIGER). Alors on peut dire, toutefois avec précau- 
tion : les ondes déterminent la probabilité de présence des particules. 

On aboutit ainsi à une conception statistique de la mécanique ondu- 
latoire, conception dont le développement formel est principalement 
dû à DIRAC, JORDAN, NEUMANN. HEISENBERG s’est efforcé d’interpré- 
ter les faits aussi bien que possible en montrant que la mesurabilité 
d’une grandeur physique a une limite absolue, donnée par la relation 
dite d’ « incertitude » 

Ap. Ag = _ 
où g est une coordonnée, p le moment conjugué, et le signe A l’in- 
tervalle dans lequel cette grandeur peut être trouvée par la mesure. Une 
mesure exacte de g exclut donc une mesure exacte de p, et, en fait, 
l’unité d'incertitude est la constante } de PLANCK. 

On comprend facilement le sens de cette inégalité dans le langage des 
ondes, Si g se trouve dans l'intervalle Ag, on doit construire un paquet 
d'ondes, qui soit considérablement différent de zéro à l'intérieur 
de l'intervalle Ag. Pour cela, on a besoin de nombreuses ondes har- 


: ART : CRT 
moniques sin —, et, puisque d’après de BROGLIE la longueur 


d'onde À est liée au moment conjugué p par la relation p — À , un 


grand Ap correspond à un grand A. BoHR a repris plus précisément 
les conséquences physiques de cette « dualité » entre ondes et par- 
ticules. 
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La mécanique des ondes ne semble donc pas conduire à des contra- 
dictions d’après cette conception statistique. 

Vous savez déjà qu'avant la mécanique ondulatoire on avait proposé 
une autre modification de la mécanique classique ; je veux parler de la 
mécanique des matrices, fondée par HEISENBERG. Il a été montré que 
celle-ci est tout à fait équivalente à la mécanique des ondes. Mais cette 
dernière s’est montrée si commode et si claire que les méthodes de 
la théorie des matrices se sont trouvées reléguées au second plan. 

La raison principale en est que la mécanique des ondes opère avec des 
fonctions continues et des équations aux dérivées partielles, quinoussont 
devenues famillières par un long usage. Mais cette similitude formelle 
avec la théorie classique ne doit pas nous empêcher de constater que les 
conceptions fondamentales de la théorie des quanta en diffèrent tout 
à fait. Les idées de Boxe sur les états stationnaires, et les transitions 
quantiques, sont caractéristiques pour la nouvelle mécanique des 
quanta. Le passage de cette théorie des éléments discontinus à la 
théorie classique des éléments continus se fait au moyen du principe 
de correspondance de BOHR, d’après lequel les lois de la mécanique 
classique sont des cas limites des lois quantiques. 

Ce fait s'exprime clairement si l’on n’emploie pas la méthode des 
ondes mais la méthode des matrices. On a observé que l’application 
des matrices n’est pas plus difficile que les opérations différentielles, sur- 
tout s’il s’agit de systèmes à spectres discontinus. En outre on 
peut de cette manière non seulement établir très clairement l’interpré- 
tationstatistique mais encore en déduire des conséquences comme la 
relation d'incertitude de HEISENBERG. 

Je veux dire quelques mots sur ce sujet. 

Nous admettons, avec HEISENBERG, que tout système mécanique 
peut être caractérisé par une matrice 


VETUICR 
GE (ane ce. | = (Qnm): 


correspondant au développement en série de FOURIER d’une coor- 
donnée g classique : 


= Da, n—vE 
k 
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A chaque indice # correspond un terme ou valeur de l'énergie W, et 
à chaque couple des termes (ou transition) #, »# une fréquence 


Lu ne 
Van = } (Wa — Wh)- 


En ajoutant aux éléments de la matrice des facteurs harmoniques par 
rapport au temps guet on voit que les éléments de la diagonale- 
Gm sont les valeurs moyennes de la quantité g dans l'état ». 

Cette simple constatation de HEISENBERG est notre point de départ 
pour l'établissement des lois statistiques. Avant de les traiter, nous 
devons dire quelques mots de la notation employée dans le calcul 
des matrices. 


Le calcul des matrices 


Nous prenons, comme définition de la somme et du produit de- 
deux matrices : 
(A no Bhux = Ann + Ban 
(ABjeu = D AuBim- 
à 
Nous notons par 2* le nombre conjugué d’un nombre complexe = 
(et non pas par :, parce qu'en physique le tiret désigne la valeur 
moyenne). 
La matrice adjointe d’une matrice À — (Aux) est 
At = (AE) = (A*.). 
On a la règle 
(AB)+ — BtAt. 
Une matrice s'appelle matrice adjointe à elle-même ou hermitique- 
quand À — A+. 
La matrice réciproque d’une matrice A est la matrice représentée: 
par A-! qui satisfait à la relation 


AA AR ST, 


où I — 


ss... 


est la matrice unité. 
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On a 
{AB} — BA", 
Une matrice U s'appelle « unitaire » (ou orthogonale hermitique) si 
DUr= Tr 
On peut aussi interpréter les matrices comme des opérations qui 


portent sur des vecteurs x ayant des composantes x,, x, (générale- 
ment complexes). Dans ce cas Ax est le vecteur de composantes 


(AX)à = > Are 
k 


Le produit de deux vecteurs x, y est défini par 


at 
(x, y) = De = (y, x)*. 
n 


La forme hermitique D'Axtx, appartenant à la matrice À peut 


nnt 


être écrite de deux manières : 
(Ax, x) = (x, ATx) ; 
car à cause de 


AGP AE 


nm 


On à 


(x, Atx) — Dr.(Dats.) = Z(Darx.) 


m 


== DAS) — DANAL 
— 22...) AY, x). 


Comme longueur d’un vecteur on définit la grandeur 
4 = V(x, %). 
On appelle l’ensemble de tous les vecteurs x qui ont une longueur finie 
(malgré le nombre infini des dimensions) un « espace d'HILBERT » ; 
dans cet espace on fixe une métrique en supposant l’invariabilité de la 


longueur. Alors toutes les transformations permises des coordonnées 
sont unitaires : car si x — Ux! est une transformation unitaire on a 


(4,24) (Ur Ur) = a UEUT) = (4; #1), 
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et en vertu de la réversibilité de la transformation, on conclut 
UTU= UUL =. 
Une forme quadratique (Ax, x) se transforme par la substitution 
x = UY! 
en 
(Ax, x) — (AUXx!, Ux!) — (UTAUX!, x!) = (UT'AUXx", x’), 
ou bien la matrice À se transforme en 
At UV" Ur 
Soit F (4, ds...) une fonction donnée, dont les variables indépendantes 
sont des matrices ; on a pour une transformation unitaire des va- 
riables ge = U-lqU 
FE = F{g'igr ee.) = UT'E(giigs : - JU. 
Les fonctions F les plus importantes sont surtout celles qui sont 
hermitiques lorsque les matrices variables le sont également ; elles 
s’appellent fonctions adjointes à elles-mêmes. 


On peut réduire toute matrice hermitique À à un système d’axes 
fondamentaux ; il existe toujours une matrice unitaire U telle que 


DA NUE (adm). 
ls Les ax sont les valeurs propres ou caractéristiques de A. Si on écrit 
AU — aU 


on peut faire opérer les deux membres de cette équation : sur un vec- 
teur arbitraire x°, et si l’on pose 


LAN C2 
on à 
Ares ar 
Les vecteurs qui satisfont à cette équation sont les vecteurs propres 
ou caractéristiques de A. 
Dans la théorie d'HEISENBERG par exemple les Wa sont les valeurs 


propres de la matrice H — (Hum), qui représente l'énergie du système 
dans un système de coordonnées q°, p°. 


Si nous représentons des grandeurs physiques par des matrices, 
il sera commode de regarder toutes les matrices transformables 
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l’une en l’autre par une transformation unitaire comme représentant 
une même grandeur. On peut appeler cet ensemble de matrices un 
tenseur de l’espace d’HILBERT ; les matrices sont les représentants 
du tenseur dans les différents systèmes de coordonnées. 

On peut aussi caractériser les points de l’espace d'HILBERT par des 
fonctions continues. Si (x1, x2, …) est un vecteur et 4(g), Pa(g), … un 
Système orthogonal complet de fonctions de la variable g, de sorte que 


J len(g) lg = à, ‘ en(g)en(g)g — 0, nn Æm 
on peut construire la fonction 
X(Q) = Xiva(qg) + Hapo(g) + ce... : 
elle représente le vecteur x dans un système de coordonnées ayant un 
ensemble continu d’axes. En effet, si l’on écrit la variable g comme 
indice on a 
n 
et cela signifie une transformation unitaire de matrice 
U — (%w), 
dont les éléments ont un indice continu et un indice discontinu. 

Cet indice g continu peut être considéré comme le numéro des va- 
leurs caractéristiques d’une grandeur physique (comme dans la théorie 
d'HEISENBERG par exémple l’indice discontinu # est le numéro des 
valeurs propres de l’énergie). Il y a un nombre infini de grandeurs 
physiques qui ont des valeurs caractéristiques 49 — w(g), fonctions con- 
tinues de g. Parmi cet ensemble de grandeurs la plus simple est celle 
qui a les valeurs caractéristiques a — q ; cette grandeur est notée sim- 
plement par le symbole g, et la même notation est aussi employée 
pour le tenseur correspondant. 

Quand on veut distinguer les valeurs propres du tenseur on écrit 
avec DIRAC pour les valeurs propres g' ou g”. On peut donner une 
représentation du tenseur g au moyen d’une matrice symbolique, 
introduite par DIrAC. Pour cela on définit la fonction symbolique (g) 
de la variable réelle g, qui est nulle pour chaque valeur de g excepté 
g = 0, où elle devient infinie mais de manière que 


2 + 


1% à(g)dg = z. 


SPA = 


M. BORN 


Le sens exact de cette définition de DIRAC est le suivant : d (g) est 
la limite d’un système de fonctions (4), (9), … tel que 


Co 
lim gag) — 0, pour g2<0, lim f f(g)e(9)d9 — f( 
n —}00 n—}00 & = 


pour toute fonction f(g). 
Maintenant nous pouvons voir qu’une représentation de la grandeur 
qg est donnée par la matrice à indices continus 


Qu = (UT —T); 
car on a 
Ja — drag — dx(a) 
ou, plus brièvement 
gx =qx. 

Cela veut dire : l'opération sur le vecteur x indiquée par g est la 
multiplication de x par la valeur propre g° de g. 

Quand g est une coordonnée d’un système mécanique on peut se 
demander quelle est la représentation du moment conjugué # dans le 
système d’axes fondamentaux des g. 

p est défini par la reletion de commutabilité 

h 


2ri 


bg — gp — 


Cette équation est satisfaite par l'opérateur 


ha. 
P— ox dq : 
car on a pour un vecteur arbitraire x (g') 
hifi, , 4 \ h 
(bg — qgb}x — 2x dg 1 AT) Lattes 


Pour représenter p sous forme d’une matrice il faut introduire la 
dérivée d'(g) de la fonction © ; alors on a 


h 
= 7 ÿ! DO SE Nr 
Para 2T1 (g q ) ; 


car on à par une intégration par parties : 


J 3ri° (g  — g'hx(g'hdg = ES 
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Les valeurs propres p’ de p sont définies par l'équation 
h dx. 


V—= = = x < 
p 271 dq PS 
c'est une équation différentielle, qui a pour solution générale 
pa 
X — 4e 


Pour définir univoquement les valeurs propres de ’ il faut se donner 
l'intervalle dans lequel g peut être choisi ; cela dépend de la significa- 
tion physique de g. Quand g’ peut prendre une valeur réelle quelconque 
il n’y a aucune restriction pour D’; mais quand get g + 274 déter- 
minent la même position du système mécanique, la fonction x doit être 
périodique par rapport à g° avec la période 27% ; on en tire 


= — - N— 0 200 2, 


? 


2e 
. ETS D. # # . 
En tout cas les fonctions e / sont les éléments de la matrice 
unitaire (à indices continus ou discontinus) U,,, qui transforme le 


système d’axes de g en système d’axes de D. 
Il y a quelque difficulté formelle à écrire la relation UU+ — 1 sous 
forme d’intégrale ; car on a 


271 


2 on 27 hog" ON 7 
furrurias ferrre Her Ar aps 
intégrale qui est divergente pour l'intervalle d'intégration 
pee, pour g—d#0. 


Mais on peut substituer à cette intégrale la valeur moyenne 


(22 x 
271 f n 
HAN HI me {1 pourg —g" 
im — — : 
Pn2Ae | o pour g Z q': 
a 
c'est-à-dire qu’on considère le système infini à spectre continu comme 
cas limite d’un système fini à spectre discontinu. 
On peut traiter le calcul des opérations différentielles comme 
ee h où 
| 271 2q 
de la même manière et avec les mêmes notations que le calcul 
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des matrices. Si A est un opérateur différentiel quelconque qui trans- 
forme une fonction x(g) dans une autre Ax(g), on peut définir le 
produit scalaire par une intégrale : 


= | ##(g)-Ax 
L'opérateur adjoint de À est défini par 
(Ax, x) ce (x, A*x), 
c'est-à-dire 


Ÿ 2*.Axdg — | x-(A+x)*dg. 
Par exemple l’adjoint de 


? = D EANE est pr = ?, 
parce que 


hCax * 2rt dx* 
PRÉDAE fx. (Prx)Fdq = fx CE dg = — fx —— 4 


et par intégration partielle 


apte) = far À ag = [2 pudg — (px, 2) 
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pour toutes les fonctions x(q) qui sont nulles à l'infini. 

Cette supposition est généralement nécessaire pour le domaine des. 
fonctions x(g) qui peuvent être considérées comme représentant des 
vecteurs dans l’espace d'HILBERT. Car pour cette classe de fonctions, 
la multiplication des opérateurs différentiels correspond d’une façon 
univoque à la multiplication des matrices à indices discontinus. 

On définit la matrice de l'opérateur 


AD, 9 = A( 2, a) 


à l’aide d’une série de fonctions orthogonales #,(q), 2(q) 


pu. 


+ 


Av,(q) — > ArnnPm(Q), 
k 


dont les coefficients ont la valeur 


Ann = [ +*-Ao,dq = (Avr, Om) ; 
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ce sont les éléments de la matrice À — (Am). Si B est un autre 
opérateur, et B— (Bmr) la matrice correspondante, le produit des 
deux matrices a pour élément 


(AB) == PAT. = D'Am(Bon, ox) = (Be, d'Asb), 
k 


LE 
k k 
puisque 
GE, )) = (0%,9) (4, C9). 
Mais 
Ja in AA 
alots 
(AB). = (By, D Atpi) — (Be, Aty.) ; 
k 


et si les fonctions +, appartiennent à la classe dont nous venons de 
parler, on en tire 
(AB }rn _— (AB, Om) 

Cela veut dire que la matrice AB correspond au produit des opéra- 
teurs AB. 

Si F(p, g) est une fonction adjointe à elle-même des variables 9, p 
d’un système mécanique, leurs fonctions propres x(g) et leurs valeurs 
propres f sont données par l’équation différentielle 

(D, g)x(g) — fx(a). 


On peut passer facilement à la représentation par matrices en posant 


x(9) = Ÿ rxen(g), 


Ad 
" 


où les #h(g) sont un système orthogonal. 
Intoduisant cette série dans l’équation, on a 


DanE(p, gJentg) — 1 Anpn(g) : 

ñn n 
si nous multiplions cette équation par #,(g) et intégrons par rap- 
port à g, nous en tirons 


Der. — Î Km 


" 
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on = JF, G)en-dg = (Epns Em). 


Cette équation est équivalente au problème de 1a réduction de la 
matrice F à des axes fondamentaux. 
On voit que les deux formes de calcul : opérateurs différentiels 
et matrices, sont tout à fait équivalentes. 
Quand F(p, g) est la fonction d'HAMILTON H{(?, g) du système mé- 
canique considéré, l'équation différentielle de SCHROEDINGER 
h à 
(a, Br 5 @ — Wx(q) 
et l'équation de HEISENBERG 


DH WE 
n 


sont identiques, c’est-à-dire sont des formes différentes de l'équation 
tensorielle 
Hx = Wx. 


Cela posé il suffit de se borner à la considération d’une des deux 
formes. Nous préférons ici le calcul des matrices parce que nous vou- 
lons démontrer que l’interprétation statistique de la mécanique quan- 
tique peut être développée d’après la théorie originale d'HEISENBERG 
d’une manière très simple et naturelle. 


1. — L'interprétation statistique de la mécanique quantique. 


Après cette digression mathématique retournons à notre problème. 
Nous considérons une fonction d'HAMILTON H(p,g) des coordonnées 
gx et des moments fx, qui obéissent aux relations de commutabilité 


h 
Digi — Qui = Zn ht) Dxbi — DID NO: 


qxdi — QiQx = 0. 

On sait que la solution des équations du mouvement 
dg_0H ‘ds = oH 
dt op, dt  oq, 


est équivalente au problème de trouver une matrice unitaire U qui 
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mette la matrice H sous forme diagonale, ou bien de trouver des coor- - 
données 


PEU U: OUT qU. 
telles que la matrice 
H(P, Q) = U 'H(?, q}U — W — (Wiômn) 
soit une matrice diagonale. 


Si F(p, q) est une fonction (adjointe à elle-même) il existe toujours 
une matrice U telle que 


CE (8), 


Imaginons que F soit l'énergie d’une partie du système donné, fi,f,,... 
étant les valeurs possibles (ou caractéristiques) de la grandeur F. 
On a alors 


PDU AE OU Un 
L 


D'après la définition donnée au début, la fvaleur moyenne dans . 
le temps, de F dans l’état », est 


ie = SU Une 
k 


La condition UU+ = 1 ou U-1 — U+ veut dire 
_ de + 
= =. 1 Le LPS 


et de Îà on tire : 


à | U 
L 


dhr Ur |? = de 
k 


2— TI 


ul 


Imaginons qu’une première mesure faite sur le système considéré, 
nous ait seulement appris qu'il est dans l’état quantique #, et 
que l’on doive déterminer par une autre mesure la valeur de F. 
La variation de F dans le temps se traduit par le fait que si l’on 
répète à plusieurs reprises la même mesure on obtient des valeurs . 
différentes f1, fa, 

Nos formules font comprendre l'interprétation suivante : Le nombre - 


Wnrx — | Urn è 
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donne la probabilité pour que la mesure nous fournisse exactement 
la Aïème valeur caractéristique de F (pendant que l'énergie du système 
total a sa nième valeur caractéristique). La valeur moyenne dans le 
temps Fun se présente alors comme valeur moyenne statistique ou 
«espérance mathématique ». À partir de là on peut aisément dévelop- 
per la théorie statistique. D'abord on peut étendre ce qui a été dit 
des grandeurs H et F à deux grandeurs quantiques quelconques A et 
B, qui sont fonctions adjointes à elles-mêmes des fx, ge; Pr, gx étant 
donnés comme matrices, A(p,q) et B(p,q) sont aussi des matrices. 
Si l’on transforme les coordonnées, c’est-à-dire si l’on prend pour 
nouvelles coordonnées 
Pr = UTPiU, gx = UT'qU, 
A et B se transforment en 
A'=U-'AU, B'=U-BU. 
I1 existe toujours un système de coordonnées Z(!) dans lequel À ala 
forme diagonale 


AU — (anônm) 
et un autre système ©(?) dans lequel B a la forme diagonale 
BE — (bô). 


Soit U la matrice unitaire qui sert de transition entre X(!) et X(à) ; 
alors où a 


A® = UAUU = ( DUstaUn), 
k # 


Bi DEN NN TURRUePT 
n 
Bien que les éléments de À, B ne correspondent plus en général à des 
fonctions simplement périodiques, nous pouvons définir B{} comme 
valeur moyenne de B quand la valeur a de À est connue. 
Alors on a 


© \ . o 
A) — DE UP, MERE 
n 


Ensuite on appellera tn — | Uinl?la probabilité de trouver la valeur 
de B, quand A a la valeur 4 ; parallèlement im — | Un |? est aussi 
la probabilité de trouver la valeur a de À, quand B a la valeur bn. 
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Quand on fixe 4 on connaît doncla répartition des systèmes d’un 
ensemble statistique parmi les valeurs 2,4, de B ; cette répartition 
est donnée par 


Wii = FU: A Wys — | Uxs Ke Rate 
Quand on connaît la kïème figne de la matrice 
CU 
U1 Us... 
De Nbre 
DU 


ou, ce qui est la même chose, la kïème colonne de la matrice 
AU Ur 
UM Ur 
QE Em , 


on connaît la répartition. 

Nous allons en donner un exemple. Soit À une coordonnée continue g, 
B le moment conjugué correspondant p. Alors —comme nouslavons 
déjà vu— ona 


La probabilité d’une valeur propre p' de p, quand g a la valeur 
donnée est donc 


Wap! —= [U, 


1p! [è —= 1, 
c’est-à-dire que chaque valeur p a la même probabilité. 

Dans le cas des variables continues on ne peut pas normaliser la 
probabilité de façon que la somme des probabilités soit l’unité ; on ne 
peut parler que des valeurs « relatives » de la probabilité. 

Cet exemple de variables conjuguées représente un cas extrême ; en 
général les probabilités ne seraient pas constantes mais fonctions de 
la valeur de la variable considérée. 

Retournons aux systèmes à spectre discontinu. 

Dans l’espace d’'HILBERT le passage d’un système de coordonnées > 
à un autre système ©’ correspond à une transformation unitaire 


AE ou AUS 
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où bien 


eee NN 
Xn = Nue == Dre 
k RME 
Quand on prend pour x le vecteur unité qui coïncide avec l'axe fon- 
damental de À appartenant à @r, 


= 000" OLO0:"-1 
(où le 1 occupe la A" place) les composantes de ce vecteur danslesys- 


tème ©’ sont 


1 * 
Xn — LVyrne 


Choisissons en particulier pour , 2’ les axes fondamentaux apparte- 
nant aux grandeurs À et B, et appelons c un vecteur fondamental 
de A (de longueur 1); on a pour ses composantes dans le système. 
fondamental de B : 


CU 


La probabilité de trouver une valeur + de B, quand la valeur carac- 
téristique de À appartenant au vecteur c a été mesurée, est 


W, = {ec |? : SaPer 
n 


On a supposé ici que À n’est pas dégénéré, c’est-à-dire que À n’a pas 
de valeurs fondamentales coïncidentes. Si on a un ensemble statis- 
tique tel que pour chaque système une certaine grandeur À prenne 
une valeur déterminée nous l’appellerons (d’après NEUMANN) un 
«cas pur » (cela veut dire le cas le plus simple possible) ; il est carac- 
térisé par le vecteur de c, et a la propriété que la statistique de toute 
grandeur B peut être construite en projetant c sur les axes fondamen- 
taux de B et en élevant au carré les composantes de c. 

La moyenne statistique d’une grandeur quelconque B peut donc 
s’écrire : 


C1, Ca. étant les composantes de c dans le système fondamental de- 
B,ona 


. 
Ÿ Brncn cn bic, 


PAUL 
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c'est-à-dire que les produits bxcz sont les composantes du vecteur Bc 
dans ce système de coordonnées. De 1à on tire 


B GT (Bec, C) = DAT 
kn 


La valeur statistique d’une grandeur B dans le cas pur c est donc 
simplement la valeur de la forme hermitique appartenant à B pour 
les composantes de c dans le système fondamental de B. 

Le cas pur représente dans la mécanique des quanta la délimitation 
la plus étroite qu’on puisse imposer à un ensemble statistique. 

Tout autre ensemble est un mélange de cas purs. Voici ce que cela 
signifie : On peut caractériser le cas pur, au lieu du vecteur c, par un 
tenseur de l’espace d'HILBERT : 

Fe cCr, Fr 260 


qui n’a jamais de valeurs négatives. On a 
B = YBuncnct = D Binlnx = SP(BF) 
kn kn 


la « trace » Sp(A) étant définie par 
SPIA) = Ÿ Ann 
ñn 


représente la somme des éléments diagonaux d’une matrice ; il vient 
à Rs Fu 
SE) = DFrn = Doonct = Ylcnl? = 1. 
ñ ñn ñn 


F s'appelle le fenseur de probabilité. 

Pour chaque ensemble qui n’est pas un cas pur, on a pareïllement un 

tenseur de probabilité (positif) ® ayant les propriétés 
SP(P) = 1 
Sp(+B) — B. 

Quand on a un certain nombre d’ensembles qui sont des cas purs, 
avec les tenseurs de probabilité F4, F®,.. on peut construire en les 
mélangeant dans les proportions w®, w),... un ensemble que nous 
appelons un mélange, ayant le tenseur de probabilité 


p — Yumi. 
n 


= ZE 
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On peut démontrer que tout ensemble peut être représenté dé 
cette manière. Enfin on a le théorème: La condition nécessaire et 
suffisante pour qu’un ensemble soit un cas pur est qu’on ne puisse le 
construire autrement que par des ensembles qui soient également des 
cas purs. 

Ces résultats sont équivalents à la restriction de la mesurabilité 
des grandeurs physiques qui dans la mécanique des ondes est 1llus- 
trée par la construction de HEISENBERG. On ne peut d'aucune ma- 
nière déterminer plus précisément un ensemble de systèmes que par 
la mesure d’une valeur caractéristique a: d’une grandeur physique À 
non dégénérée. Toute grandeur B non commutable avec À (p. ex. 
coordonnée et moment conjugué) n’est jamais mesurable en même 
temps que A. Car cette impossibilité d’intervertir le produit veut dire 
que le vecteur c appartenant à 4x ne coïncide pas avec un axe fondamen- 
tal de B; donc les composantes c;, ©, de c dans le système fonda- 
mental de B ont des valeurs qui sont certainement différentes de 1000..., 
01000...,… Ia probabilité wx — | cx |? d'existence d’une valeur caracté- 
ristique x de B est donc déterminée et différente de o et de 1. Doncil 
est inconcevable qu’une valeur caractéristique déterminée # de B et 
aucune autre ait toujours été trouvée en même temps qu'a. 

D'ailleurs on peut aussi très facilement déduire la loi de réciprocité 
d’'HEISENBERG pour la restriction de la mesurabilité des grandeurs non 
commutables. De l'égalité dérivée ci-dessus (Ax, x) — (x, A+x) on 
déduit pour la valeur moyenne statistique de AA+ dans le cas pur c : 


AA — (AAte, c) — (Atc, Ate) > 0, 
parce que le carré de la longueur d’un vecteur n’est jamais négatif. 
Soient A,B des grandeurs hermitiques, À — A+, B — B+ et À un 
nombre réel ; on a 
(À + àB) (À + AB)t — (À + AB) (À — ÀB) > 0, 
Ai + LB? — {AB — BA) > 0 


pour toute valeur réelle de À. Par conséquent on a 


A-B+ AB BA Zo. 


On remplace maintenant 
À par ôAæA—AÀ, 


B par ,iB = 8H 5: 
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alors Comme on peut le voir facilement, il vient 


dA-0B — BA — AB — BA 
donc 


SAE.6B: > — 43E pa 


En particulier, quand p, q sont des variables conjuguées, 


h 
pq — gp 7 
on a 
et DE 
NI ne MEL k 
FA OT 


Si l’on définit maintenant 
, Ap = V5’, Ag = Vie, 
on a en définitive 
ap À, 
et on reconnaît tout de suite que c’est la limitation la plus exacte pos- 
sible de la grandeur d'incertitude. 

Cet exemple montre nettement que la méthode des matrices où 
plus exactement des tenseurs de l’espace d'HILBERT permet d'exposer 
la mécanique quantique d’une manière simple et claire. Les opérations 
différentielles et intégrales de la mécanique des ondes ne sont impor- 
tantes que lorsqu'il s’agit des phénomènes ondulatoires propres, 
c'est-à-dire de systèmes à spectres continus dans lesquels les parti- 
cules viennent de l'infini ou disparaissent vers l'infini. J'en donnerai 
des exemples dans les conférences suivantes. Mais même dans ce cas 
il est souvent plus commode de rendre d’abord le spectre des valeurs 
propres discontinu par un changement des conditions et de le trans- 
former à la fin en un spectre continu par un procédé de passage à la limite. 

En toute rigueur, dans la mécanique quantique le temps est aussi 
un opérateur (tenseur de l’espace d’'HILBERT) comme toute autre 
grandeur physique et en effet elle est conjuguée à l'énergie H — E du 
système. 

11 y a donc une restriction pour la mesure simultanée de l'énergie 
et du temps ; 


0 
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La détermination d’un état stationnaire (avec une énergie exacte) 
demande donc nécessairement un temps long ; cela n'a donc aucun 
sens précis de parler des états d’un système à un moment déterminé. 
Mais comme approximation il est souvent avantageux de faire 
abstraction de ce couplage entre l'énergie et le temps. On a alors le 
droit d'employer le temps £ comme paramètre. Au lieu de l'équation 


U-'H{p, q)}U — W (matrice diagonale) 
ou 
HU — UW—=o 
on a l'équation différentielle 


et ceci doit être interprété ainsi : On introduit dans H{(#, g) pour fe, ge 
les matrices #x(0), gk(0), qui correspondent à la valeur £ = © ; puis 
on cherche une solution U(é) de l'équation qui tend pour £ = © vers la 
matrice unité. Si l'on pose 

Eat) — UT'(AR(O)U(E), Pat) = UT(#BA(O)U(#), 
on obtient les matrices des coordonnées et des moments conjugués pour 


le moment t. Toute grandeur A(f, g) est maintenant aussi fonction det 
et l'ona 


A(t) — A(p(t), g(t)) = U{(#A(#(0), g(o)}U(é) — U-(#Ato)Uf#). 


Dans la notation employée jusqu'ici A(o) et A(é) sont des grandeurs 
physiques différentes ; car elles correspondent à des tenseurs ayant 
des axes fondamentaux de directions généralement différentes. Mais 
il s'agit d'autre part de la même grandeur physique parce que les lon- 
gueurs des axes fondamentaux de A(o) et A{(é) sont égales, donc la 
réduction aux axes fondamentaux de A(o) et de A(é) donne la même 
matrice diagonale : 

A0) = AXE) = (tndum). 
Dans le système appartenant à A(0) on a donc 
At) = U-1(#A%o)U(f) = U-(AAKAU(), 
et de là on déduit 
A?) = U(8 AU (à. 
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Ici on a relié unitairement deux représentations d’une même gran- 
deur A(f) (« méme » dans le premier sens). Donc on peut employer les 
conceptions statistiques développées auparavant et on trouve le ré- 
sultat : 


Vu = | Us (f) | = | U:,(5) |? 


est la probabilité de trouver exactement la valeur 4x pour la grandeur 
À au temps { si À avait eu la valeur 4, au temps { — 0. 

On parle plus particulièrement de probabilité de transition dans le cas 
suivant. Soit un système mécanique d’énergie H,(p, q), influencé par 
une force perturbatrice extérieure, de manière que l'énergie totale soit 
donnée par 

H(D, 9, t) — H,(9, q) + (9, q, t). 
Pour 2 grandeur À on prend l’énergie H, du système non perturbé. 
On introduit dans celie-ci1 les solutions 2°, g° du système non perturbé, 


c'est-à-dire les matrices pour lesquelles H,(p°, 9?) — W' est une ma- 
trice diagonale ; et on détermine la matrice U par l'équation 


_ + Hf°,g,U—0,  U(o) = 71. 


Alors 

27 — | Uz, (6) [E 
est 12 probabilité de transition n — k pendant le temps é. On peut aussi 
remplacer l'équation différentielle pour 1a matrice U par un système 
d'équations différentielles ordinaires en choisissant un vecteur cons- 
tant arbitraire # et en introduisant le vecteur x — Ux? comme indé- 
terminée. On a pour ce dernier vecteur 


h di. RAI 2 
271 dt Date. Lt 
ou plus explicitement 


h dr, 


& 
LE D + Din = 0. 
mm 


Nous rencontrerons cette dernière forme d’équation dans les appli- 
cations. 
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2, — La théorie quantique de la désintégration radioactive des noyaux 
atomiques. Remarque sur la catalyse chimique par absorption. 


Une des plus belles applications des méthodes de la mécanique 
quantique est la théorie de la désintégration radioactive des noyaux, 
donnée à la fois par M. Gamow et par MM. CoNDON et GURNEY. 

La théorie classique rencontre une grande difficulté à combiner les 
faits mis en lumière par les expériences de RUTHERFORD sur la disper- 
sion des particules z traversant des atomes lourds avec les résultats 
des mesures de la vitesse de particules 4 émises par les même atomes. 

De la répartition des particules 4 qui rencontrent un atome d’Ura- 
nium on peut déduire que la loi de COULOMB est valable jusqu'à une 
distance de 3,2.10=!% cm. du noyau ; l'énergie électrique du système 
noyau-particule z à cette distance a la valeur : 


(e—2:4,77-107, E — 92-4,77-107%). 


Mais l'énergie de la particule 7 émise par l'Uranium lui-même n'est 
que 6,5.10—$ erg. 

En imaginant un modèle du noyau, nous pouvons construire une 
courbe pour le potentiel de la force entre le noyau et la particule x. 
Cette courbe a la forme d’une hyperbole jusqu’à la distance 3,2.10-2 cm ; 
ce point dépassé elle s’incurve et, très près du novau, elle tombe brus- 
quement de telle manière que la particule z puisse osciller autour d’une 
position d'équilibre stable. L'énergie de cette vibration est la même que 
l'énergie du rayon % après l'émission, représentée dans notre figure 
par la ligne droite horizontale. 

Comment la particule z peut-elle s'échapper de cette région inté- 
rieure du noyau et traverser la montagne de l'énergie potentielle ? La 
théorie classique ne donne pas d'explication de ce phénomène. 

Mais la mécanique quantique peut surmonter cette difficulté et en 
même temps donner la déduction d’un autre fait frappant de la 
radioactivité, la loi de GEIGER et NurraLL. Cette loi établit une con- 
nexion entre l'énergie FE des particules 7 émises par un noyau radio- 
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actif et la constante À de désintégration (c'est-à-dire la probabilité 
d'émission par seconde) ; elle a la forme 


log À — 4 + bE. 


Les valeurs de l'énergie des particules x diffèrent les unes des autres 
pour les différentes substances d’une série de désintégration au moins 
dans la proportion de 1 à 2 tandis que les constantes de désintégration 
diffèrent dans la proportion de 1 à 10%; ce fait est rendu par la 
formule précédente. 

La mécanique des quanta explique le fait de la désintégration 
radioactive et en même temps la loi de GEIGER-NUTTALL. 

Ilest essentiel de remarquer que cette explication est fondée sur le 
fait qu’une barrière d'énergie n’est pas un obstacle absolu au passage des 
ondes. Cette circonstance est bien connue en optique. Lorsqu'un rayon de 
lumière tombe sous l'incidence rasante sur la surface d’une lame à faces 
parallèles, il est réfléchi totalement, comme onle sait. Sil’on approche une 
autre lame de verre de la première à une distance qui est comparable 
à la longueur des ondes lumineuses, un peu de lumière entre dans 
cette deuxième lame. C’est que, même dans le cas d’une réflexion 
totale, il y a un mouvement d’onde dans le second milieu. C’est une 
onde qui se propage le long de la surface et dont l’amplitude diminue 
proportionnellement à une fonction exponentielle de la distance 
normale à la surface. Après quelques longueurs d'ondes elle est deve- 
nue imperceptible. 

L/onde ne transporte pas d'énergie de la surface à l’intérieur du 
second milieu. Mais si la seconde lame de verre est rapprochée de la 
première à tel point que la vibration dans l’air la rencontre encore 
d’une manière perceptible, la réflexion cesse d’être totale, et un peu de 
lumière pénètre dans la seconde lame. 

Les ondes de DE BROGLIE se comportent de la même manière dans des 
régions de haute énergie potentielle. Elles franchissent un tel obstacle et 
cela d'autant plus facilement que la fréquence y est plus grande, c’est-à- 
dire que l'énergie Av d’une particule est plus grande. Cette idée a déjà 
servi auparavant à expliquer divers faits comme par exemple l’ex- 
traction d'électrons des métaux par des champs électriques très 
intense à des températures relativement basses. 

Gamow, et, après lui, quelques autres savants ont calculé l'émission 
de rayons % en étudiant une onde diminuant avec le temps, onde qui 


» 
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se trouve à l’intérieur du noyau, et qui, pénétrant dans la barrière- 
d'énergie, se propage après sa sortie comme une onde sphérique à 
l'infini. 

Cela nous amène à employer des opérateurs non hermitiques avec 
des valeurs caractéristiques complexes. 

Les conceptions de la mécanique des quanta expliquées dans les. 
premières conférences ne peuvent être appliquées sans avoir au préa- 
lable posé des hypothèses spéciales. C’est pourquoi il me semble dési- 
rable de développer cette théorie avec les conceptions éprouvées dela mé-. 
canique des quanta. C’est de cela que je vais brièvement vous parler. 


Nous considérons le spectre des valeurs propres d’une fonction po-- 
tentielle à symétrie sphérique V(r) de la forme d’un cratère, repré- 
sentée dans notre figure. Pour éviter les difficultés des spectres continus, . 
imaginons le noyau entouré d’une très grande sphère qui limite le 
champ Ÿ extérieurement ; alors la fonction d'HAMILTON H — + V 
(T' est l’énergie cinétique) n’a que des valeurs caractéristiques discon- 
tinues. Parmi celles-ci, celles qui sont plus petites que la hauteur du 
bord du cratère, c’est-à-dire que le maximum V, de V, se divisent vi- 
siblement en deux classes : d’une part celles qui sont sensiblement- 
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écartées les unes des autres, et qui sont relatives aux fonctions propres 
oscillant dans l’intérieur du cratère et tendant rapidement vers zéro 
à l'extérieur : ce sont là les valeurs caractéristiques et les fonctions 
caractéristiques «intérieures ». Il y a d'autre part celles qui sont 
extraordinairement rapprochées et se rapportent aux fonctions propres 
qui oscillent hors du cratère et qui, à l’intérieur descendent rapidement 
à zéro : ce sont là les valeurs propres « extérieures », et les fonctions 
propres « extérieures ». 

Admettons donc maintenant que V(r) se compose avec une 
grande approximation d’un potentiel attractif et d’un potentiel ré- 
pulsif. Soient V.(r) et V,(r) des fonctions analytiques définies pour 
toutes les valeurs de 7, de telle sorte que 


VAT NV ANNE DOTE AE TS; Vin  Z 0, 
Vi(r) >> V:(r) » PE ry, Vz'(r) < 0. 


Soit ensuite le véritable potentiel 


Vi(ro) = Va(ro) —= Vo, 


Va POUPEE, 


Ne, Va(r) » ET: 


Nous allons maintenant considérer deux systèmes auxiliaires avec 
les fonctions d'énergie 


H=T+V;, HT: 
Admettons que leurs valeurs et fonctions propres soient : 


Win Um. pour Hi, W,,7%, pour Ep, 


de sorte que nous avons l'identité : 


(H + Wan )Um — O, (Æ + Wu 6; 


Les valeurs propres de H,, où H, concorderont approximativement 
avec les valeurs propres intérieures, où extérieures de H. H, et H, 
ne sont pas commutables, ni par conséquent simultanément mesu- 
rables. Si une mesure de H, donne la valeur propre Wn, cela signifie 
que la véritable fonction d'énergie H a une valeur proche de W» et 
une fonction propre qui partout concorde presque avec #n ; c’est-à-dire 
que la particule x est à l’intérieur. Alors elle ne peut être en même temps à 
l'extérieur, c’est-à-dire que la solution correspondante de l'équation 
d'onde (H, — Wjv — o est v — 0. De même la mesure d’une valeur 
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-W, de H,, signifie que la particule z est dans l’espace extérieur, donc 
que (4H, — W)x = 0 doit être satisfait par # = 0. 

Si l’on a de cette manière établi au temps £ — o que la TE d. 
est à l’intérieur (Wa, #» pour H,, v — o pour H), il peut arriver qu’on 
trouve au temps f la particule dans l’espace extérieur (W,, v, pour H,, 
# — 0 pour H,) ; et à la vérité cela arrivera d'autant plus fréquemment 
que la sphère extérieure limitative sera grande, c’est-à-dire que le spectre 
de Wy sera serré. Tant que la sphère est bornée, on obtient une 
probabilité périodique de transition D»,(#), c’est-à-dire que la parti- 
cule x oscille entre l’espace intérieur et l’espace extérieur (nous 
-considérerons un tel cas plus tard en parlant de la catalyse de 
contact chimique). Avec un espace extérieur infiniment grand dy 
s’annule toujours, sauf quand W — W, et dans ce cas le rapport 


æ 
w = =. tend vers une limite déterminée quand # augmente indéfini- 
ment, c’est-à-dire qu’il existe une probabilité de transition par unité 
de temps. 


La miseen œuvre mathématique de cette méthode se fait comme suit. 
L'équation de SCHRÔDINGER pour le phénomène véritable est : 


(# — e )vE == 10} 


Nous cherchons une approximation de la solution par une somme : 


2 274 ÿ 
Vi) —Ù an(te ' Un FD b,(t)e EU 7 


m F 


dans laquelle ne doivent se trouver que des valeurs d'énergie Wn, 
W,. qui sont petites en comparaison de la hauteur du cratère V,. Si l'on 
introduit cette expression dans l'équation différentielle, qu'on mul- 
tiplie celle-ci par #,* ou par v,* et que l’on intègre par rapport à 
l’espace, on obtient en vertu des relations d’orthogonalité et de 
normalisation de #n, Vy. : 


d e UA QE 2TV 1l D be 
(D 


h 
2ri re 
a 


= 
ok nl Vu (a) TV nl 25 2TIV PR 0 
(à +25 mm © mx g }+De Ci sa b He ; 
nr 
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Ici sont utilisées les abréviations suivantes. D’abord on a introduit les 
fréquences 


va = ÿ (Wn — Wa), “ti 


Puis on forme des fonctions 


FO (») = V(r) TE. V,(r) pour #7 F0) 


(o) » PET 
a SEE { 0 pour #<#, 
EU l AGE V,(r) » RUE 


d’où l’on tire les éléments de matrice (dS — dxdyd2) : 


[F0 uvas =, [Frs = Fi 


LA 


{a) + (a) Te pe 
fr UnutdS = Fi EF unotdS = 5% 
enfin nous posons : 
* OR 
| vu) 19 = mure 


Ces grandeurs sont sans exception petites, mais d'ordres diffé- 


rents. Les plus grandes sont encore relativement FO et EG) ; dans 


Fr par exemple, le produit FO ; qui dans l’intérieur est sensi- 
blement différent de zéro, est multiplié par la fonction v, qui à 
l'extérieur oscille et à l’intérieur est presque égale à zéro. Mais 


dans FQ et FC) les différences des énergies potentielles sont, dans les 


régions où elles diffèrent de zéro, multipliées chaque fois par deux fonc- 
tions propres, qui dans cette région sont très petites (elles disparaissent 
comme des exponentielles) ; nous négligerons donc ces termes. De 
même nous pouvons omettre les termes contenant des :, parce qu'ils 
ne donneraient que des corrections d'ordre supérieur. Alors on obtient des 
équations de l’espèce introduite d’une façon générale à la fin du $ 2. 
Les coefficients am, b, forment les composantes d’un vecteur de pro- 
babilité, leurs carrés donnent les probabilités de transition. 
Nous partons maintenant de l’état initial 


am(O) __ LEAAS b,, (0) AS 


L2 
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qui signifie que, jusqu’au temps { — © la particule a été observée à 
l’intérieur avec l'énergie War. Ensuite nous introduisons ces valeurs 
initiales dans les termes dont nous venons de reconnaître la petitesse 
(termes du premier ordre) et tirons de la seconde équation : 


ne (a) omiv 1 
— == — e 
5 bd. = F 1€ c 


L'intégrale qui s’annule pour { — o est : 


y 
CLOUS AE 
D pe er 


n nÀ 


La probabilité de transition d’un état intérieur # à un état exté- 
rieur À est donc : 
À : DT à Sin? v,l 
Pa | be (FR) ER one 
ce qui est en fait une fonction périodique du temps. 
La probabilité que la particule saute de l’état intérieur à un état exté- 
rieur quelconque est la somme 


LA n PRES 


Si nous passons, en augmentant indéfiniment le rayon de la sphère, 
au cas limite de valeurs propres extérieures WA réparties d’une façon 
continue, la somme se transforme en une intégrale : 


2 
L sin + (W, — W;}é 
= (ryef apr [Ra 


oo 
: SE 


Soit maintenant À — À(E), une fonction qui découle, à la limite, de 
l'échelle à degrés équidistants À par interversion de l'équation 
W = E ; soit \'(E) sa dérivée. Nous posons : 


et 
x =} (E — Wah; 
alors on a : 
2x\? ht . hx sin x\? 
“= (7) =] eme) ) 
— 75 Wn 
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donc : 


: AU 47° 
NE HE W, > 
{—> 0 l } À ) 
C’est 1à la probabilité de transition par unité de temps, c’est-à-dire 
la constante de la désintégration radioäctive. Le calcul de cette cons- 
tante se réduit à la détermination de l'élément de matrice : 


- 20 


12% = | Ft ds fa | Vilr) — Vi(r) É uVIAS. 


On peut encore simplifier cette expression en utilisant la symétrie 
sphérique de la fonction potentielle. Alors #», v, deviennent des pro- 
duits de fonctions du rayon 7, Rn(r), Pu(r), par des fonctions sphériques. 
Si nous nous restreignons à l’état de moindre énergie de la particule 
liée au noyau # == 0, la fonction sphérique correspondant à R, sera 
constante et égale à r ; des relations orthogonales des fonctions sphé- 
riques il résulte que : 


a 


FO | (Vi Vi)RiPurdr. 
L “o 
J'ai effectué le calcul explicite de cette intégrale pour deux hypo- 
thèses différentes sur la fonction potentielle. Dans la première le po- 
tentiel à l’intérieur est supposé constant avec une montée brusque: 
vers la valeur du potentiel de COULOMB, pour 7 = 7, (voir la figure) ; 
la fonction V, est alors ainsi définie : 


= ï pour 7 Fo, 
Vo pour 7 . Vo; 
tandis que V,(r) est égale au potentiel de COULOMB : 


: 2e?Z 
VA ILSe Re 

La seconde hypothèse aboutit au modèle suivant : au centre se trouve: 
une charge négative — Zi ; celle-ci est entourée vers l'extérieur par 
une charge superficielle positive, plus grande, Zas de sorte que l'effet 
d'ensemble vers l'extérieur répond à la charge positive : Ze — (Za=— Zi)e 
(voir la figure 4). 


Les fonctions propres intérieures peuvent être données exactement. 
dans les deux cas. 


\ 
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Pour le calcul des fonctions propres extérieures j'ai utilisé un procédé: 
qui a été indiqué indépendamment par BRILLOUIN et par WENITZEL : 
on introduit dans l'équation d’onde pour # —7 P(r) 


Ê Br? : 
+ TUE — Vu = 0 


l'expression 


Fig. 4. 


et on prend la fonction y comme inconnue ; alors on obtient : 


LEE se 
re À == 2M(E Va) ya 


PES 
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En négligeant le terme multiplié par k on obtient : 
y = + V2M(E — V:) 
d’où l’on tire : 


2Ti 5 
2% | VY2ME — Vijdr 
P(r) = 4 € : 1 


RS 
nm: 

Cette formule représente suffisamment le caractère de la solution, 
-surtout la chute exponentielle vers l’intérieur. Pour la normalisation il 


faut y ajouter le facteur _, dans lequel I, est le rayon de la grande 


VL” 
: sphère ; celui-ci finit par disparaître du résultat, parce que les valeurs 
propres de l'extérieur se comportent asymptotiquement comme 
h? 


EM" 


.d'où l’on tire 


SC at HU ENEe 
Les Va ue 
Les intégrales qui apparaissent dans le calcul ont la forme : 


Ro 1 f(r)e-4\ dr 
oÀ a ER à ? 
-dans laquelle f(7) est une fonction rationnelle, g(r) la fonction transcen- 
.dante qui figure dans l’exposant de P(r). On la traite le plus simple- 
ment d’après la méthode suivante, souvent employée : 

Une contribution essentielle à l'intégrale n’est fournie que par les 
alentours de la position 7, où g(r) a un minimum ; on développe donc 
g(r) d’après des puissances de 7 — 7 et on ne retient que les termes 
jusqu’au second ordre ; alors l'intégrale peut s’évaluer. 

Le résultat dans le premier cas (potentiel constant à l’intérieur), 
peut se ramener à la forme suivante : 


(4 2 = 
af IVÊÈTE = E )dr 
w p(E, Vo) € sons 


“ici g(E, ro) signifie une fonction lentement variable de l'énergie E de 
la particule 4, et qui ne dépend, en dehors de constantes (charge élé- 
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mentaire «, masse M de la particule 7, numéro .atomique Z du noyau) 


2 
que du rayon 7,; 7, est la distance pour laquelle Ni, 


Cette formule diffère à la vérité de celle de GHIGER-NUTTALL, mais 
fournit une relation de même allure que cette dernière entre w et E. Elle 
est aussi très semblable aux expressions données par GAmow, CONDON, 
GURNEY et autres savants (!). Un de mes élèves, M. ANDRESS a examiné 


numériquement cette formule. L'intégrale de l’exposant se laisse facile- 
ment évaluer et conduit à la fonction : 


arc cos = — Her MP 
2) 

ANDRESS a donc examiné quelles valeurs il faut admettre pour », 
pour que dans une série de désintégrations (série Ur-Ra) les w et E ob- 
servés résultent de la formule. Pour cela il faut encore avec ATKINSON et 
HOUTERMANS faire de petites corrections aux valeurs d’énergie obser- 
vées, surtout à cause du choc en retour que le noyau émetteur subit. 
On trouve avec une erreur probable d'environ 1 °/, : 


/ 


Va a 


LT Ur II 10 Ra RaC' RaA 
010!? cm Lt? TPE 1,05 1,0 1,10 0,98 


J1 en résulte donc que tous les rayons sont presque constants et 
égaux à 10- cm. Ainsi l’idée fondamentale de la théorie est justi- 
fiée. 

On peut sans doute appliquer la méthode exposée ici à beaucoup 
d’autres problèmes encore. Je voudrais en citer un: la catalyse chimique 
par adsorption dans des corps solides. Au dernier congrès annuel de la so- 
ciété de BUNSEN, dans lequel la question de la catalyse forma avec 
d’autres le centre des discussions, HABER a émis l’hypotèse que des 
phénomènes de la mécanique quantique (par exemple échange d’éner- 
gie par résonance) y jouent un rôle essentiel. Mon ami FRANCK et 
moi avons essayé de nous faire une représentation définie du phéno- 
mène. Je veux dire ici quelques mots de ce travail. 

La catalyse sur les corps solides se produit, autant que nous le sa- 
chions, à la condition que les deux partenaires de la réaction soient ad- 
sorbés. En outre elle ne s'effectue manifestement pas sur toute la sur- 


(x) R. d’E. ATKINSON et F.HOUTERMANS, Z5. f. Phys., 84, 656, 1929 ; J. KuDAR, Zs. f. Phys., 
53; 166, 1929. 
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face, mais seulement en certaines régions « actives » qui sont probable- 
ment caractérisées par des inégalités de structure ou par la présence 
d’atomes étrangers dans le réseau. Enfin elle semble s'effectuer d'autant 
mieux que la distance des atomes de la substance du catalyseur est plus 
petite Il s’agit de comprendre pourquoi, précisément dans de telles 
conditions, une énorme accélération de la réaction se produit. 

La lenteur de réaction des gaz repose sur l'existence reconnue par 
ARRHENIUS d’une «chaleur d’activation ». On peut aujourd’hui d’après 
LonpoN et HEITLER la comprendre au point de vue de la mécanique 
quantique. 

Quand par exemple deux molécules biatomiques s’entrechoquent, 
l'énergie du système (en négligeant l'énergie cinétique des noyaux) 
peut être considéré comme fonction de leur distance 7 ; cette fonction 
a alors d’après les calculs de LONDON l'allure suivante (voir la 
figure) : 


V 


Ice: 


Il existe un minimum, qui correspond à la position d'équilibre #, 
du produit de réaction ; mais avant celui-ci se trouve un maximum ”;, 
et qui empêche les molécules, quand l’énergie cinétique est faible, 
puissent se réunir. Seul le très petit nombre de molécules dont l’énergie 
relative est plus grande que le maximum près de 7,, peuvent parvenir 
à la réaction. C’est pourquoi celle-ci s'effectue d’une manière très lente 
dans l’état gazeux. 

Or, d’après les conceptions courantes, l'effet du catalyseur doit être 
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de diminuer la liaison des molécules qui doivent être amenées à la 
réaction, ou bien de les dissocier et d’adsorber les atomes isolés : de cette 
manière le maximum (qui est précisément l'énergie d'activation), doit 
être abaïissé ou supprimé. Reste à savoir si ces suppositions sont vrai- 
ment propres à expliquer le processus catalytique. En tous cas nous 
croyons qu’en outre un effet spécifiquement quantique s’y ajoute qui 
repose sur ce que précisément dans la mécanique quantique une bar- 
rière d'énergie n’est pas un empêchement absolu pour le mouvement de 
particules. A la vérité pour l’état gazeux il n’existe aucune probabilité 
pour qu'un passage ait lieu, parce que les temps des chocs sont beau- 
coup trop courts. 

Cette probabilité est une fonction approximativement périodique 
(au moins immédiatement après le commencement de l'interaction) de 
la forme : w — sin? 27vé. Quand le temps de l'interaction est court, re- 


lativement à la période T — , on a: w — (2rvf)? quand il excède 


cette valeur T — £, on trouve w d’ordre 1. Nous voyons dès lors l’effet 


de l’adsorption précisément en ce qu’elle prolonge ces temps d’interac- 
tion. Les deux partenaires de la réaction séjournent pendant longtemps 
à la surface, manifestement, d’une manière particulièrement fermeaux 
centres actifs ou centres de désorganisation, et ils sont d'autant plus 
serrés que la constante du réseau de la base est plus petite. Alors ils 
ont le temps de passer dans une nouvelle position d’équilibre d'énergie 
plus faible, même quand une barrière d’énergie considérable est pla- 
cée devant. A l’état gazeux cette probabilité reste petite même quand 
on augmente la pression, c’est-à-dire quand on augmente le nombre des 
chocs ; car pour chaque couple d’atomes la probabilité part dela va- 
leur zéro et augmente d’après la loi : w — (2rw)?, et le résultat par ad- 
dition des probabilités, pour un grand nombre de couples d’atomes, 
n'est pas aussi grand que celui qui résulte d’une augmentation du 
temps f. Pour évaluer le temps T, nous avons pris comme modèle le plus 
simple une courbe d'énergie avec deux minima, dont l’un avec 7 — 0 et 
V — o répond à l’état avant la réaction, et l’autre avec 7 — a et 
V — — U répond à l’état après la réaction (voir la figure 6). 

Pour plus de simplicité nous composerons cette courbe de deux 
paraboles. Alors nous pourrons appliquer nos considérations avec de 
faibles modifications. Seulement nous ne considérons pas le processus 
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dans le temps, mais l’oscillation stationnaire qui satisfait à l'équation 
aux valeurs propres : 
(H — E)Y = 0. 
Nous écrivons la solution comme auparavant sous la forme : 


(7) =Ÿ anum(r) + D buvulr). 


ni nr 


V 


Fig. 6. 


Alors on obtient pour les constantes a», b, les équations 


Le Np EU) — 
alW, — Ë) APE, — 0. 
U 


GMA Dante = 0) 

d 
(dans lesquelles nous écrivons maintenant les indices 1 et 2 au lieu de 
1 et a). Nous admettons maintenant qu'il existe deux niveaux d'énergie 
égaux de l’oscillation non couplée. Ce n'est naturellement jamais 
exactement le cas pour les molécules isolées ; mais justement à cause 
de l’adsorption sur la base, chaque état d’une molécule sera considé- 
rablement indéterminé de sorte que notre supposition est admissible. 
Soit donc 
Wi= W, = Wo 


& 
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L'oscillation couplée aura dans la première approximation (« pertur- 
bation séculaire ») la forme : 


BY by — 0, FR 
et une énergie : 
E == W, ia É: 
où EF’ est petit. Puis on obtient : 
XÉ' — FC y — — 10; 


JE! — FSx — 0, 
d’où l’on tire la séparation des termes d’énergie : 
BE —(/riel. 
A celle-ci correspond une fréquence de battement : 
A = } (Ex —E); 


donc le temps du battement + par rapport au temps moyen d’oscilla- 
tion des systèmes non couplés : 
LT SN ho 
Too Er ne 
y L Ar FU QE > 


a 


Or nous avons évalué les éléments de matrice pour les paraboles 
d'énergie : 
Vire 2m ur? Va = 27°vp(r — a) — U; 


ici ,, », sont les fréquences des oscillations non perturbées. Nous intro- 
duisons des grandeurs sans dimension en posant 


y Vo 
V— vi + va, a=., P=, (a + P= 1), 
U V 7 
P —\ 47" ua)", u = p’ o = re Ps = A 
L 2 
où L. est la masse «effective » de la particule vibrante ne ne , 
2 


et nous supposons que 


5 KI, U «KI. 


M. BORN 


Alors 


Les fonctions propres de l’oscillateur sont 


ne ie T 4T?uv 
9 — #94 | ue — 
Un N H,(é)e 2 ; Ny \/: u | v ds À = ne 


où HA(£) est un polynôme d'HERMITE ; pour les vibrations à gauche 

(potentiel V,) et à droite respectivemen(V,) on a au lieu de » et de £ 
Pa 

Vis ë, = V Tes et Va 


Alors on a 


(1) Pa fre l_/Pa P£ ) Le; 
1554 e N,Ny En | pr e) H, (W (@—1))F(o)l de, 


ou 


et FC) est la même intégrale entre les limites p, et co. 


Si on écrit la fonction à intégrer sous la forme e— tr ce ne sont 
que les valeurs de ® près du minimum qui contribuent à la valeur de 
l'intégrale. On peut considérer comme constants les deux polynômes 
d'HERMITÉ dans cette région, parce que le point 2% est hors de l’inter- 
valle où les polynômes oscillent. Nousles approximerons d’une manière 


très grossière par les puissances les plus hautes, Hn © (2Ë)*. 
Nous posons 


où p est une valeur moyenne proprement choisie. 
Alors on a 


"00 _Wy 
Pon = C [7 pe FT à 
—@ 


On ne peut pas remplacer F(p) par une valeur moyenne parce que 
F(+) s’annule pour o = Po (donc F(2) — ). 
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On définit maintenant la fonction 
J2 
ge) = 7 fe) — Log F(+) 


et on cherche les minima de (s) près de p, ; il y en a deux, l’un fa 
à gauche, l'autre De à droite de Po et 


g(u)=0,  v(p)=0 
Alors le développement de la fonction g(e) au voisinage de p, est 
pe) = (pi) + = #1): (e (4 He 


et on a approximativement 


Ne Re : — 
EU HQE P)do — Ce ?(P1) ‘ha RENE UE Nes 
CEE mmail °°) En ° 71 
et de même 
EF) — Co—v1(0:) (= 
min Do 
Enfin on trouve 
an e. V2r. El Co 
h ÿo7e Il un 


Pour l’évaluer, le mieux est de choisir ÿ comme racine de l'équation 
f'(p) = 0, c’est-à-dire 


ph RE 


L’intersection des paraboles, F(4) = 0, est approximativement 


pp VER + 7 pe — (VE he) 
pp (VIP — = po + SEE + T a), 
où 
pe ee 
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De cela on tire avec la même approximation 


h 
ei + a = RÉ 108 D ; pt, 


PA 16 
MURS 7 
Heu (x) (25 s) 
Nous trouvons le résultat : 


Su 5 
se Lo” dr a ane I 
SN a RER Te 
y. nimi 4, r16Pa ‘h 

V 25 + pp 
Û 


Quand l’oscillation à gauche à wn quantum d'énergie, m — 1, nous 


trouvons pour la proportion de la durée de battement 7 — T à la durée 


Av 
, ne Te CLR 
d’une vibration propreT — , — Æ 
4p , 60 

ee CHAN PPa (à nr. IÈaf 
a —— 2 e 
4ù Av DRE ui AE 12 

nur DAS 


Ici on peut introduire au lieu de P la hauteur de l'intersection des 
paraboles : 


2; 
V,= PT sp TP page, 
à 2e 2 FR D 
: © 
Si nous choisissons comme unité de longueur 1A — 10-° cm., 


comme unité de fréquence 10% sec., comme unité d'énergie 1 grande 
calorie par mol, nous avons 


= 0,94», — 66 pa, 


A MiM» . : F + SALRS . 
où pu — nn est le poids atomique « effectif », c’est-à-dire, dans le 


cas de deux atomes égaux la moitié du poids atomique. 

Maintenant nous nous sommes fait une représentation du phénomène 
si important au point de vue technique de la synthèse de l'ammoniaque 
NH, en nous rattachant à l'opinion que le processus primaire con- 
siste en la formation d’hydrazine 


H 
APR" à 


ÉLUS 
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L/'état initial est formé par des molécules au repos les unes à côté 
des autres 
H H 
| N=N | 
H H 
Le processus est essentiellement la rupture de deux liaisons N — N 
au cours de laquelle les”atomes N s’éloignent les uns des autres. Et ici 
c’est précisémentÆla chaleur d'activation qui doit être surmontée, 1a- 


ro F 
by 


œ 
o GS; 40 A 


Fig. 7. 


quelle dans notre représentation est figurée par le maximum d’éner- 
gie à l'intersection des paraboles. 

Nous avons maintenant à l’aide de données optiques et thermiques 
évalué grossièrement les nombres d’oscillations et la différence d’éner- 


gie U et nous avons abouti à admettre que 
Y=— 0,9; Va =="; donc = 1102, 
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tandis que U sera situé entre 20 et 60 grandes calories ; à ce dernier 
nombre correspondent des valeurs de 
PR 
— 
comprises entre 7 et 20. On obtient alors 
pour #20 HUE log — 31,64 — 9 log'° a — 8,5, 


» UE :60 07 — 20% logi° — 31,6 a? — 22 log'° a — 8,0. 


Si l’on dessine ces fonctions (fig. 7), on voit qu’elles décroissent en une- 


© 
pente très rapide et que peu à peu vers a — 0,8A elles conduisent à l’ordre 
de grandeur 14 ; 1à nous avons donc + — T-10“c 1 sec. On obtient 
donc ce résultat que, avec des durées d’adsorption d’une seconde des- 


déplacements des positions d’équilibre de 0,8 À sont possibles. Mais. 
c’est là ce qu’on peut attendre à peu près dans la formation de l’hydra- 
zine comme grandeur du déplacement des atomes N. Nous croyons, 
que ce résultat parle pour la justesse de notre hypothèse que l’effet ca- 
talytique de surfaces repose principalement sur la fixation de molécules. 
dans des positions voisines pendant des temps prolongés. 


8. — Sur la largeur naturelle des raies spectrales 


Dans la théorie classique de la lumière, on explique la largeur d’une 
raie d'émission par l’amortissement exercé sur l’oscillateur rayon- 
nant par son propre rayonnement. D’après la théorie de HER/IZ 
un oscillateur harmonique de fréquence », perd en un temps Af, grand 


< : Éd > * 

par rapport à | , mais petit vis-à-vis de la seconde, l'énergie 
0 

(1) AE = — 5 v vo P°A4 


où P représente le moment électrique du dipôle et où la barre indique 
la moyenne par oscillation. D'autre part, l'énergie de l'oscillateur lui- 
même (masse #», charge e) est 


_—. 
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par suite on a 


__ am?e?v? 
3cm 


dE, 
dt RE X E, J( 
avec la solution 
É— He 2: 


le moment P sera par conséquent représenté par la fonction 


—" no —"1 TA _ Ti\ \ 
P(é) = Pie T! cos 2rvt — 7 Pué CPR RER ST 


2 


Si on développe ceci en une intégrale de FOURIER, on obtient 


A+ 00 


Ph= |]  pte""a, 
Se où 
é Ti : Po ; e 
AA — 2m 7 LE «: 
PL) = f DU = Y — 27%) — Ÿ) y + 2700 + v) 


Dans cette formule on peut négliger le second terme, parce que», + y 
n'est pas nul dans l'intervalle o = y  « et y est très petit en com- 
paraison de »,. 

Maintenant la moyenne dans le temps de P? pendant le temps Aë 
est, entre £ et & + Aï : 


»- . nt+ At I {+ P+o 
Pr=;, | P(#)dt =; | p(y)P(t)e"""dvat. 
Î l «/ — 


L'énergie émise pendant le temps Af est donnée par (x). 
Nous obtenons alors l’énergie totale I émise en sommant AE pour 
tous les intervalles Af de £ = 0 à é = © : 


16m (Tr rivl 16n°v,* f° rivtP? 
1 je | je POP at | pra Rs. 
L'énergie émise dans l'intervalle de fréquence dy est donc : 
PERTE item 
TO Lane gr — 


I(») a son maximum aux environs de » — y, et la valeur du maximum 


est I(,)—= = Les valeurs de » pour lesquelles la valeur de I(v) est la 


moitié du maximum sont 


& MN 


On get première pour meer de le langer matarelle Be gene À. 

L'enceene Deore quantique de Foie net. por expliquer ke Rae 
ere, lier eos ke comte de froquence de FER, où RË 
RER en Ve RR Dre RER. GRR LRQ BETERRE RÉ 
RRQ SERRE FL ET ER RER RETIRE RRQ LE ie RIRE RENTRER 
Les ncatats qu'en ebteet par là ROUE es QRRRR BR PRRVERE 
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QUELQUES PROBLÈMES DE MÉCANIQUE QUANTIQUE 


On obtient la répartition de l'intensité de la manière suivante : Si 
la raie est une raie de résonance, son niveau final — l’état fondamental 
— doit être considéré comme rigoureusement défini. L/indétermination 
de la fréquence ne peut être produite que par les fluctuations du niveau 
initial. L/intensité I(v)dy dans l'intervalle de fréquence dy est donc- 
analogue à (2) et donnée par : 

Y dy Ch 
RQ) + r: 


(3) I{v)dy ==> 1 


1, est l'intensité totale de la raie. Ce résultat est d'accord avec l’ex- 
pression classique. 

Soit une raie produite par le passage du niveau À au niveau B ayant 
les valeurs respectives »4 et 5. I/intensité relative dans l’inter- 
I(»)d» 

I 


— de cette raie, est égale à la probabilité. 
0 


de trouver y4 —- /B entre y et y L dy. 

S1 les valeurs 4 et »B oscillent autour de leurs valeurs moyennes 
5 et % d’après la loi (2), la probabilité de trouver le niveau A 
entre vâ et y + dy£ et la différence v4 vB entre y et y + dy est. 
donnée par 

j ELA dys IE a dr : 
Tr 4r?(va — V64)t + ya nm 4m — y — VB)? HE y? 

La probabilité cherchée résulte alors de l'intégration par rapport: 

à toutes les valeurs de »4 et nous obtenons : 


valle de fréquence d : 


a al LE cs Es 
I(v) Ve d : = [4r?(va — voÂ) + 4°] [ar?(va — y — v,B)? Et Y5°] 
(4) 
| = J, Her 
A NN ER) 

On parvient à ce résultat par une intégration complexe autour du. 
demi-plan supérieur de »4. y, — vi — % est la fréquence à laquelle on. 
s'attend d’après la condition de fréquence de BOHR (). 

La raie À —> B a donc la même forme que la raie de résonance. Sa 
largeur est la somme des largeurs du niveau initial et du niveau final. 
Nous montrerons maintenant qu’on obtient les mêmes résultats par la. 
mécanique des quanta. 


(1) Il est permis, d'intégrer de — à +, car les y sont très petits en comparaison: 
de vhet À. 
; 0 Ô 
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Les équations de MAXWELL pour le vide sont : 


rot + LE 2 0, div H= 0 
otH— 1 ; divE=— 0 


On introduit un potentiel scalaire d et un potentiel vecteur À : 


il est clair qu’on peut toujours choisir A, de façon que ® — 0. Alors 
on a 


et on peut encore soumettre À à la condition que 
. TE 2 
div A = 0 pour {= 0; 


Les équations de MAXWELL signifient maintenant 


À 
—+ = + 
AA — 5 = 0, div A = 0, 
où 
= 5 a+ y? à + =. 


est l'opérateur de LAPLACE. L//énergie contenue dans le volume V est 
20 = fc + H?) dx dy dz 


rer on me 
= & | nr + (rot A)° À dx dy dz. 


e 

Considérons maintenant les oscillations propres d’une enceinte vide 
de volume V. 

Comme conditions aux limites, nous exprimons que les a 
tangentielles de É et les composantes normales de H sont nulles. Cela 

+ 

signifie que les composantes tangentielles de A sont nulles. 

Nous considérons d’abord les oscillations propres de l’enceinte, 
c’est-à-dire nous posons 


> 


À — À: cos (2rvÉ + vs) 
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QUELQUES PROBLÈMES DE MÉCANIQUE QUANTIQUE 
et les équations canoniques du mouvement : 


dp, _.  oK dq, __dH, 


FAR dt op, 
Rappelons en outre la formule asymptotique bien connue donnant 
le nombre des oscillations propres dans un intervalle de fréquences dy : 
4V#d» 


4 
Nous allons maintenant transformer ces formules classiques en for- 
mules quantiques. 
Nous conservons l'expression de la fonction hamiltonienne mais nous 


y remplaçons les ,, gr par des matrices qui satisfont aux règles de 
commutabilité canoniques 


he 
Pr — hi or Pbr — Phi, Gas — 94 = 0. 


De plus, nous conservons les équations qui relient les g et les À ; 
dans celles-ci, les composantes de A% restent des fonctions ordinaires, 
les composantes de À deviennent au contraire des matrices ; par suite 
les composantes de É et de H deviennent également des matrices ; la 


dérivation par rapport au temps qui intervient dans le calcul de E est 
donnée de la manière bien connue par 


YA PES Mr, 
SE = + (HA — AK). 

Cette transposition a naturellement la signification d’une hypothèse 
physique, tout à fait analogue à celle que l’on fait quand on transpose 
par exemple l’expression classique de l’énergie d’un atome d’hydro- 
gène dans la mécanique quantique. 

Les valeurs fondamentales et les fonctions fondamentales d’un os- 
cillateur sont connues. On peut en déduire d’une manière simple celle 
d'une enceinte vide. Pour le s-ième oscillateur, on a les valeurs 


fondamentales : 
I) 
W,= him + ;), Ms — 0,1,2, - -- 
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et les matrices des coordonnées (en utilisant l'expression de la 


masse Li) : 


More e. à De K 
/ M? ; (e +1) 
Qliss Ps + 1) V pa (Hs + I) e\# 


qe + 1, %) —\/ M à (ns + I)e 


qg{nns) = 0 lorsque #:—“# EI; 


l'énergie totale de l'enceinte devient donc : 
10 — Ÿ W.; Ca Wa + Ÿ hvsns 
ee 2 


S S 


+ IN ; - à ; é 
où Wo — Dh représente une énergie au zéro absolu. Celle-ci est 
S 


infinie, ce qui est sans doute une difficulté de principe pour la théorie ; 
mais elle ne nous gênera pas dans les considérations suivantes qui ne 
concernent que les différences d'énergie à partir de l'énergie du zéro 
absolu. . 

Les matrices gs du système d’oscillateurs équivalent au champ 
de rayonnement, s’obtiennent à partir des matrices des oscillateurs 
individuels par une opération que j'appelle la « fusion »: 


\ 


Qu, Na, cs M, M 0) 


SR NS à [AN = 
— Ôns, n° * Ôno, ne! On, nr Gtts, Ms )0n. ru, M'osrs 5 


leurs éléments sont donc toujours nuls, sauf lorsque tous les #,, #2... 
restent inchangés à l'exception du seul n, qui varie de + 1. 

Supposons maintenant qu’un atome de fonction hamiltonienne x 
se trouve dans le champ du rayonnement ; supposons le noyau de 
l'atome infiniment lourd et fixé au point x,, y, &. Admettons enfin 
que seules jouent un rôle les longueurs d'onde du rayonnement qui 
sont grandes par rapport aux distances moyennes des électrons et du 
noyau dans l'atome ; nous pouvons alors prendre pour calculer l’in- 
teraction de l'atome et du champ de rayonnement les valeurs des inten- 
sités des champs au point x,, Yo, 4, Où se trouve le noyau. 

Si ds représente l'énergie du Hohlraum seul, &,! celle de l’atome 


ns. = 


dèrie un db ‘à ré in 
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isolé, nous obtenous pour le système total enceinte vide plus atome, 
y compris leur interaction, une fonction hamiltonienne 


BH = Hg + Hs + U 


P'représente le moment électrique de l'atome. 

Nous traiterons ce système par les méthodes du calcul des perturba- 
tions. Comme système non perturbé, nous prendrons la fusion des deux 
systèmes : atome et champ de rayonnement. La fonction perturba- 
trice U est une matrice dont les éléments dépendent des indices n 
des valeurs fondamentales de l'atome et des indices n,, n2,.. ms... des 
valeurs fondamentales du champ ; on obtient 


Un; m,hn ns, nn ;nm...,n,---) 
TTC MT 
=; WAm,n,, M, :--,Mm,mMm; M, :--). 


Dans cette formule, nous pouvons laisser de côté la sommation par 
rapport à s, puisque pour chaque système d'indices n,, m2... ; nj/, 2. 
11 y a au plus un seul des g, qui soit différent de zéro. 

A partir de maintenant, nous désignerons en abrégé par s le système 
d'indices n,, M, m.. et nous écrirons pour l'élément de la matrice 
perturbatrice : 

PS - 

Nous pouvons supposer que les éléments diagonaux Py de la 
matrice P sont nuls ; cela signifie que l’atome n’a pas de moment élec- 
trique moyen. D'après cela les propriétés essentielles de la matrice 
perturbatrice U sont les suivantes : U}*, n’est différent de o que 
inzn et si tous les m — n pour toutes les valeurs de s, sauf une 
seule pour laquelle m = mm +1. 

Les équations démontrées à la fin du $ 2 et qui donnent les proba- 
bilités de passage, sont : 


ns 
CEST 114$ = O. 
271 di PC. 
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Soit maintenant Frs une valeur fondamentale du système non 
perturbé 3, + 4, ; nous posons : 


Ans) ta tje 


et obtenons après une transformation simple : 


LS N LU En Enoy 
gi On T Zn É ARE 

Telles sont les équations sur lesquelles DIRAC a édifié sa théorie 
des probabilités d'émission et d'absorption. Ainsi l'émission corres- 
pond à la solution pour laquelle au temps { — 0, tous les 4 sont nuls 
sauf un seul, 4 — I, qui se rapporte à l’état » de l'atome, avec ab- 
sence de rayonnement (s — 0, c’est-à-dire M, — M9 = … hs = .… = O). 
Pour trouver cette solution, DIrAC utilise un procédé d’approxi- 
mation assez primitif : Il substitue dans la somme les valeurs initiales 
données des ans et peut alors effectuer facilement l'intégration. Mais 
cette approximation n’est valable que pourunintervalle de temps très 
court après l'excitation, d’une manière précise, pour un intervalle 
petit vis-à-vis de la vie moyenne de l’état excité de l'atome. Après un 
temps aussi court, l'énergie émise est, d’après la relation d’indétermi- 
nation de HEISENBERG, tellement indéterminée que la variation d’in- 
tensité à l’intérieur de la raie spectrale ne peut pas être observée. Si 
on extrapole le résultat de DIRAC pour des temps très longs, on obtient 
toujours, même en tenant compte les approximations plus élevées, des 
raies spectrales infiniment fines. Cela n’a naturellement aucun sens 
physique, car dans ce cas, les approximations ont depuis longtemps 
cessé d’être valables. La largeur des lignes ne peut donc être obtenue 
que par une solution qui soit valable avec la même approximation 
pour toute la durée du rayonnement. 

Une telle solution est encore intéressante à un autre point de vue : 
elle donne les probabilités des niveaux atomiques individuels en 
fonction du temps, tandis que les calculs de DIRAC ne déterminent ici 
aussi que le premier moment. C’est ici que commencent les calculs de 
MM. WEIssKOPF et WIGNER qui donnent une solution plus rigoureuse. 
Un résultat important de leur étude est, que les passages individuels 
sont à une bonne approximation indépendants les uns des autres, 
comme on l’admettait dans l'ancienne théorie quantique. Les proba- 
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bilités satisfont, lors de l'émission spontanée, à l'équation de désin- 
tégration radioactive 


de #2 >à Dita — > Wyfnis 
: ; 


où les ‘x représentent les probabilités de passage du -ème au k-ème 
état. 

Nous considérons le cas particulièrement simple d’une raie de réso- 
nance, c’est-à-dire une raie pour laquelle il n’existe qu’un seul passage 
spontané partant du niveau B, l’état final étant l’état fondamental A. 
Si l’on tient compte des propriétés des éléments de matrice ue ’être 


toujours nuls sauf lorsque n >< n’, et qu’une des différences ns — ns’ 
est justement égale à + 1, on peut ranger les 4ns d’après les passages 
possibles de la manière suivante : 

En premier lieu, on a à considérer les grandeurs 4% qui sont initia- 
lement égales à 1 (atome éxcité, pas de rayonnement) ; ensuite, les 
grandeurs 4,, où dans s un seul des indices #, est égal à 1, tous les 
autres étant nuls (atome dans l’état fondamental, un quantum de 
lumière dans le rayonnement). Ensuite viendraient les combinaisons 
d'indices qui se déduisent de cette dernière par le saut + 1 d’un indice; 
ces combinaisons ramènent, soit à 450, Soit à des états dans lesquels 
le rayonnement est caractérisé par un système d'indices s composé de 
zéros et de deux unités (atome excité, deux quanta de lumière dans 
le rayonnement). 

L'hypothèse principale de WEISSKOPF et WIGNER qui nous permet 
de résoudre les équations est que l’on peut négliger ces états et 
tous les processus qui s’y rapportent. Car ces états seraient produits 
par le processus suivant : l’atome passe de l’état inférieur 41, à un 
autre état excité, tandis qu'un quantum de lumière est encore émis 
(et non pas par exemple un quantum de lumière absorbé). Bien que 
la théorie en elle-même autorise de tels processus, on sera cependant 
en droit, pour des raisons physiques, de les considérer comme telle- 
ment improbables qu’on puisse les négliger dans les équations. Natu- 
rellement il faut alors montrer que la solution ainsi obtenue remplit 
les conditions avec une grande approximation. 

Nous négligerons encore les états 4,,, dans lesquels la fréquence du 
quantum émis diffère, strictement parlant, de celle qu’on attend par la 


= — 
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condition de Boxer. Cette omission est aussi légitime que l'autre. 
Elle n’est pas absolument nécessaire pour la résolution des équations 
mais elle est plus commode pour le calcul (?). Il nous reste donc main- 
tenant les équations suivantes : 


h daso ; N'y, % Œpo — East 
Re | = 0, 
2ri dit PE ” 
2" = \ 
EE LU NOT 


\H—0,S = it 
NS — à 
ne doit être prise que pour les À qui correspondent à un petit intervalle 
de fréquence autour de », (la fréquence de résonance), de y, —: 
à , + «.« doit être petit devant »,, mais grand en comparaison de la 
largeur de la raie. 

Nous écrivons pour simplifier : 


L'indice À signifie la combinaison d'indices la somme 


430 = b, ar = 4, UÙUy = = Un 


et introduisons, au lieu des énergies, les fréquences : 
Ero — Eur = A(vo — vx) 


où », est la fréquence de résonance ; x est alors la fréquence de la 
kme oscillation propre. Nous obtenons alors : 


le RUE 
D +ŸUe tv vol 
k 


Æ + 27 — V0) 
FE LE U; e 27 (Vi — Volt b. 
Nous posons la solution : 


vi 
b—=e : 


(1) Si l’on ne fait pas cette dernière omission, non seulement on obtient la largeur des raies 
mais de plus le milieu de la raie est déplacé de la position prévue d'après Bohr (v. Dirae, L. e.). On 
ne peut calculer ce déplacement qu’en connaissant exactement la fonction perturbatrice U pour 
toutes les fréquences. Cependant il faut observer que notre U n'est valable que pour les longueurs 
d'ondes grandes par rapport aux dimensions atomiques. Approximativement le déplacement 


est grand en comparaison de la largeur de la raie, mais cependant petit en comparaison de 
la fréquence moyenne vo. 
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et obtenons : 


pe REZ VE — 
ia = U,e [Y + 2m — V)l. 


? 


en intégrant 


Uz lente voR nd 1]. 


y —= : : 
LY — 2TU(Vz — Vo) 


en tenant compte de la condition initiale 4x — 0 pour £{ — 0. 
En substituant dans la première équation différentielle, on obtient : 


2e | U, 1e lex nee On 0 on A 


À — 27% — vo) 


Si on fait tendre l'enceinte vers l'infini, la somme par rapport à k se 
transforme en une intégrale par rapport à la variable » — w. | Ur |? est 
une fonction continue de y, de sorte que nous pouvons définir une 
fonction (>) par : 


> RU = NAi0, d= y, = Y} 
ee À Le 


Il vient alors : 


Det et il Her 10) reel, 


Nous pouvons facilement calculer l’intégrale d’une manière appro- 
ximative. Comme d’après l'expérience y est très petit vis à vis de w, 
on peut sans grande erreur considérer, dans l'intégrale, /(v) comme 
constant et égal à f(v,), car l'intégrale possède pour y — », un maxi- 
mum très aigu. 

Il est permis d'intégrer de — © à ne o , car choisissant € 5 y les 
parties de — co à—+setde+esà+os rdlent 

Nous obtenons alors 


Vo+E f(v) —ifiy or) | 7 Ro 1h PETER rt 
AE a lv=2iç) Le dif (ve) 
et par suite 

Y = (vo). 


La valeur que nous avons posée pour b est donc bien une solution : 
l’état excité diminue progressivement suivant une loi exponentielle. 


M. BORN 


On peut calculer la valeur f(v,) de la fonction pertubatrice, que 
nous avons posée comme DIRAC : 
BO_ 2x1 here 
Le FUar — … Te 1 PÇAÏ)g — ee * nt P Aï)g, 


En y substituant l'expression de gk(nx, nr + 1) pour #%x — 0 et en 
tenant compte de la valeur donnée plus haut pour le nombre des 
oscillations propres dans l'intervalle dy, on obtient : 

4: V 
or 


fopde = TE | PA [# d 


Le ER EE 


En prenant la moyenne par rapport à toutes les positions de l’atome 
dans l'enceinte et en utilisant les propriétés d’orthogonalité et de 


+ 
normalisation des Af, on obtient : 


2rv, 2 ES 
10 = Sp Pol 
La probabilité de passage par unité de temps est donc : 


CAR : Gary 
À — ( Ê 5 1e = 2} = 2r{(%) = RE | P, |? 


ce qui coïncide exactement avec le coefficient de probabilité A 
d'FINSTEIN. Si l’on prend pour | P, |? la valeur qui correspond à l’éner- 
Rues 

— Brmv; 


e , . . Es 
gie au zéro absolu d’un oscillateur harmonique, | P, |? on 


8 2,,2 
A2Ÿ— +. ce qui n'est pas autre chose que la valeur classique de 
rappelée au début. 
Pour la probabilité de trouver un quantum de lumière de fréquence » 
dans l'enceinte, nous obtenons : 


RO 
AT (VE — Vo) Uug 1 
pour des temps très longs cette expression se réduit à : 


: [Ur 
Ar | = 
LR are 


— 2Yl 


Ü — 2 COS 2T(vw — m}tre 


+ |: 


Fan 


Pour obtenir la probabilité par intervalle de fréquence dy, sommons 
depuis &, jusqu’à k, sur un étroit intervalle de fréquence dy. Nous 


2 
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pouvons considérer cette somme comme la probabilité W{(v)dy pour que, 
longtemps après l'excitation, un quantum de lumière de fréquence y 
se trouve dans l’enceinte ; d’après la définition de /(v) elle vaut 

Y D | 

FA — VV) + 


W(v)dv — 


Dans cette formule, on a remplacé (y) par sa valeur /(v,) au milieu 
de la raie, ce qui est permis, puisque nous considérons (y) comme 
constant dans le petit domaine de fréquence dans lequel l'intensité 
de la raie spectrale est notablement différente de zéro. Ce résultat est 
d'accord aussi bien avec l’expression classique qu’avec (3). 


= -p7--- AU 


_—__——— — 


À 


Fig. 8. 


Le cas général où l’on a à considérer plusieurs niveaux d’énergie de 
l’atome n'offre maintenant aucune difficulté de principe. Il suffira 
d'étudier le cas de trois niveaux d’énergie. Supposons qu'au temps 
{ — 0, l'atome soit excité au niveau € ; de là il peut passer en B et en 
A ; soit À l’état fondamental, que l’atome peut également atteindre 
directement par le saut € + A. 

Nous devons nous attendre à ce que la raie B > À ait une largeur qui 
ne dépende que de sa propre probabilité de passage ; au contraire la 
largeur C -»> A dépendra des probabilités de passage C -> A et B > A, 
la largeur € -> B de toutes les trois. 

Nous posons de nouveau pour abréger : 

do = C; Aak — Ak 
Apx = D; dar,t — Ai 
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où les indices À et Z, ainsi que kl, signifient respectivement que dans la 
suite s—%,,n,.. tous les », sont nuls à l'exception de 57 — I ouwm—T, 


‘ — yl . 
où #x et #1 — I. Alors nous poserons de nouveau c — € Ÿ tandis que 
b se composera de deux termes ; d’une part l’état B augmentera pro- 


portionnellement à € d'autre part il s’amortira lui-même pro- 


portionnellement à e—"*, y, y, sont des constantes à déterminer. 

En faisant les mêmes hypothèses que ci-dessus, on simplifie le sys- 
tème d'équations de manière à pouvoir le résoudre commodément ; 
mais nous ne ferons pas les calculs, nous nous contenterons d’écrire 
les résultats. Posons : 


27 327 Ve __ 32rvs? 


EVE 3h | Po (£, Y2 ES | Poa ke NN 3hc? Pr Le 


alors ‘y, est l’une des constantes cherchées et 
rtf 


l’autre. 21, 2%, 2, sont les trois probabilités de passage en concor- 
dance avec la théorie de DIRAC. 

De plus on trouve pour la probabilité de rencontrer wn seul quantum 
de lumière de la fréquence y après un temps assez long 


ETS œ 
W:(v)d = R(yat valet AR(v vs)? 


cette formule donne la répartition de l’intensité dans la raie C —+ A. 
A côté de cela, il y a une probabilité de rencontrer deux quanta de 
lumière de fréquences y et v': 


€ dvd! 
W{s, v'}dvdv = TP . 
CR TR [(ya + ve) + 4m v + vi — V0) ] [y + 47 (v — vs)?] 


Pour obtenir la répartition de l'intensité dans la raie de résonance 
B + À, il faut intégrer cette expression par rapport à y’ ; on obtient 
ainsi 


Ê dv 
W d\ ET Vis — - 5 _— 
a(v) À Y1 ie Y2 y” < 4ar°(v — V3)” 


Cette formule ne diffère pas de celle que nous avions obtenue aupa- 
ravant en ne tenant compte que de deux niveaux. 

On obtient la forme de la raie € — B par intégration par rapport 
à y. On est amené à une intégrale de la même forme que (2) ce qui 
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Wat .UThEn, LE 
Sn 7 Cr ant AT) 


C’est justement la même répartition de l'intensité que (4), équation 
qui résultait de la combinaison des deux termes flous C et B. 

La largeur de Ia raie B — À est donc 7, ; celle de C—B, y = 1 + y 
et celle de € — B est la somme de ces deux : y + 3. 

Ia méthode de calcul qui vient d’être exposée peut s'étendre sans 
difficulté au cas le plus général où un état quelconque est excité au 
temps £ — 0. Les cas déjà calculés laissent d’ailleurs prévoir la forme 
du résultat général. 


Note ajoutée à la correction 


Paragraphe 2. — Le calcul des probabilités de transition dans le cas de la catalyse par 
-contact a été amélioré récemment en tenant compte du couplage entre la molécule et le cristal 
{M. BorN et V. WEIsskopF, Zettschr. f. Phys. Chemie. Abt B., 1931). 

Paragraphe 3. — La théorie de la largeur des raies spectrales que nous avons exposée a 
été publiée par V. Weïssxopr et E. WIGNER dans le Zestschr. f, Phys. 63, p. 54, 1930. De 
plus, V. Weisskopf a traité par la même méthode le problème de la fluorescence de résonance 
(Gôttinger Dissertation, Ann. d. Phys. 1931). 
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1. — Dans tout problème de Physique mathématique, l'intégration 
des équations aux dérivées partielles qui définissent l’état interne 
présente des difficultés très inégales suivant la forme des frontières 
et les conditions imposées le long de ces frontières. 

Les propriétés générales sont assez simples pour l'équation de 
LAPLACE, et plus généralement pour les équations du type elliptique ; 
on sait en outre traiter toutes sortes de problèmes à frontières sphé- 
riques, en hydrodynamique, en élasticité, en électromagnétisme, soit 
à l’aide des diverses fonctions de GREEN, connues, soit par des dévelop- 
pements en fonctions sphériques, très étudiés, soit même à l’aide 
d'équations intégrales. Ce qui facilite le travail, c’est l’accord des sy- 
métries de l’équation de LAPLACE d’une part, de la forme sphérique de 
la frontière, d’autre part. Cet accord esten partie conservé quand on 
passe aux frontières ellipsoïdes, mais la perte de l’isotropie spatiale de 
la frontière crée des difficultés souvent sérieuses. Pour d’autres 
formes, le travail est beaucoup moins avancé et souvent très diffi- 
cile. 

On peut faire des remarques analogues pour les équations aux dé- 
rivées partielles de type hyperbolique rencontrées dès le milieu du 
xvirre siècle par les géomètres qui se sont attaqués à la théorie de la 
propagation du son (Note I).— Dans ces questions acoustiques les plus 
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simples, les conditions aux limites, en temps et en espace, se présentent 
aussi avec une symétrie analogue à celle de l’équation aux dérivées par- 
tielles On peut alors obtenir des solutions satisfaisantes soit sous forme 
finie, soit à l’aide d’intégrales, soit par des développements en sérics, 
qui sont bien connus. 

Mais les difficultés apparaissent, et considérables, dès qu’on associe 
une frontière sphérique à une équation aux dérivées partielles hyper- 
bolique, comme cela se présente dans un problème hydrodynamique 
très simple, qui va servir de thème à ces conférences. 


2. — Un globe sphérique (le globe terrestre simplifié), animé d’un 
mouvement de rotation uniforme (vitesse w), est entièrement recou- 
vert par un liquide pesant ; on suppose que la surface libre du liquide 
est une sphère de rayon R, lorsqu'il tourne avec le globe. Sous l’action 
d'un champ de gravitation externe, de potentiel P‘/”, le liquide peut 
prendre un mouvement d’oscillation contraint, avec soulèvement pé- 
riodique (27 : 6). (Note II). 

On veut étudier ce mouvement lorsque la période est plus longue 
que la demi-période de rotation : 0 < 20. 

Réduit à sa forme la plus simple, ce problème de marées se ramène 
à l'intégration de l’Equation aux dérivés partielles hyperbolique 


(A) 


oP  »P eu EN See 
CT CE ) de 


avec des conditions frontières 


X(18P, — 26P)) + y(20P; + 1:0P,') LEA LAS OUE 
$ R(4° — W) A États 


ib(P — P°) - 


sur la sphère extérieure 


+ += R! 
et 


A[50P; — 2wP,] + Y[2wP," + 10P,/]  £zP! 
R(4w? — 0?) Ne 


sur le fond sphérique 


h constant. 
Ce problème, je ne sais pas le résoudre. Je vais étudier avec quelque 
minutie les difficultés auxquelles on se heurte quand on l’attaque par 
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les différents procédés les plus naturels ; cela servira d'exemple, et 
montrera quelques-unes, au moins, des questions qu’il faudrait avoir 
d’abord bien étudiées, pour attaquer utilement ce problème, — et les 
problèmes réels plus difficiles — avec chances de succès. 


: : : Ps c : { 
Si la vitesse de rotation w était nulle ou simplement plus petite que : , 


l'équation serait elliptique, et, on sait que la solution serait unique 
et bien déterminée — c’est dans ce cas que s’est placé POINCARÉ dans 
son célèbre mémoire des Acta Mat. (t. 7), de 1885 «sur l'équilibre d’une 
masse fluide animée d’un mouvement de rotation. » 


8. — Mais, dans le cas que je veux examiner, l'équation étant hyper- 
bolique, il n’y a pas de résultats mathématiques généraux, noüs appre- 
nant si les conditions frontières déterminent une solution, en déter- 
minent une seule, ou sont incompatibles. 

Précisons cette difficulté en examinant le cas où le potentiel pertur- 
bateur a la période 2r : À en longitude : 


P,=P,ere, 


Cherchons Px en coordonnées cylindriques &«, z, 5 — Va? + y? 
Avec res variables, l’équation devient : 


d°P; 4w? — 0? °P, ns I oP, k? 
dm”? 2 d2° D 0 D? 


BP 0 


Elle admet pour caractéristiques les droites : 


2 
ann 
2 
BG AR SE 
en posant 
vaut — 6? 


(dr 


ë, n étant des paramètres variables. Notons en passant que l’équa-. 
tion obtenue en prenant Ë, 7 pour variables, soit : 


est, sous le nom d’ «Equation d'EULER et de Poisson », la mieux étudiée 
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des équations hyperboliques à deux variables (?), sans que nous puis- 
sions ici en profiter. 

Supposons pour simplifier qu’il s'agisse, dans le plan méridien, 
d'un problème à frontière circulaire unique. Dans le cas d’une équa- 


dPz 


tion elliptique, une relation linéaire entre Pr et la dérivée normale 


en chaque point de la frontière détermine en général une solution 


a sur un ar 
arc 


ouvert quelconque détermine, si elle existe, la solution dans tout le 
domaine où elle est analytique. 

Dans le cas qui nous occupe, menons les caractéristiques tangentes 
au cercle, et celles qui passent par les points de contact (fig. 1). Sur 


unique à l’intérieur. Par ailleurs, la donnée de Px et 


Het oP4 , F SEAT 
d'arc Ty, T;, on peut se donner Px et FT É On détermine ainsi une solu- 


tion unique (?) dans le losange quadrillé A,T.CQT4. Sur les arcs T\T, 
el TT, T4, on ne peut se donner que Px. En effet Px est déjà connu sur 
les segments de caractéristiques T, C, et T,C,, et ce sont là des don- 


(1) CE DARBOUX, Leçons sur la théorie des Surfaces, t. II, L. III. 
(2) Cf. E; PrcARD, loc, cit, XIIe leçon et suivante. 
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nées suffisantes pour obtenir la solution dans les losanges hachurés 
PAT, TaCaTs A4. Enfin sur l'arc T, Ts, il ne reste plus aucune quan- 
hté arbitraire disponible. La solution est en effet déterminée dans les 
parallélogrammes C,T,B:Cs, C1CB3T4 puisqu'elle est connue sur C,T, 
et CoCs, CaCs et CTs. Enfin jelle est déterminée dans le losange C;B, 
A,B; puisqu'elle est connue sur C,B, et CB. 

Ainsi, pour les seules conditions à la frontière qui aient fait l’objet 
d’une étude approfondie des mathématiciens et qu’on sache aujour- 
d’hui utiliser, interviennent nécessairement les points de contact des 
caractéristiques, dans le plan méridien, qui correspondent dans le 
problème à trois dimensions, à certains parrallèles critiques. 

Laissons donc de côté ces remarques — plus ou moins inspirées par 
les problèmes acoustiques (z remplacé par {) et cherchons directement 
une solution au problème posé, au moyen de développements en séries de 
solutions simples. Cette recherche présente des points intéressants, 
qu'il importerait d'approfondir et de généraliser. 


4. — En coordonnées sphériques, l'équation (A) est fort compliquée, 
mais permet de reconnaître aisément l'existence de solutions simples 
de la forme : 

17 Are M, "(L). 


A est une constante imaginaire, k un entier ; n pourrait être pris quel- 
conque, mais nous ne considérerons que des valeurs entières. On trouve 
pour MA une équation différentielle (Note III) analogue à celle des fonc- 
tions tessérales (harmoniques sphériques), mais beaucoup plus com- 
pliquée. 

Pour en voir la signification, prenons pour coordonnées 

10z 22 
PT pre 


LV OR == 


On ramène ainsi l'équation (A) à l'équation de L'APLACE. À toute 
intégrale de l'équation de LAPLACE correspondra donc une intégrale 
de l’équation précédente. En particulier, posons : 


2 


: à : £ 
CRE les he nm 1 EE me ER 


nous aurons des solutions de la forme : 


De 
R?T! 


AëkignTi(v) et: Be** Toi: 
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A et B étant des constantes complexes. Ces solutions s'expriment 
aisément en fonction de 7 et /1. On a en effet : 


rep de DE 
TR 4w° NE V4w? np 4w*u? l 


> N° 2 20 ni 
de 2e î = r\/z — RL V = SA TE 
BR et y sont réels dans la bande équatoriale comprise entre les paral- 
lèles + 1, imaginaires purs dans les régions polaires correspondantes. 
Pour 4 = lu, & est nul, y est infini. Ces parallèles correspondent, dans 
chaque plan méridien, aux points d’intersection du cercle avec les 
caractéristiques passant par l’origine. Ils ne s’étaient donc pas in- 
troduits jusqu'ici. 

Les solutions simples obtenues s’écrivent ainsi : 


; ENS 0 \ 
Aegika yn Le _ UE ( rar H 
(3) is ( ire Zu Var 


Ne 2 


; I A à 40 

(4) A En Sent 

I1 semble que nous ayons obtenu ce que nous cherchions. Les solu- 
tions (3) sont celles qui conviennent pour un problème intérieur, les 
solutions (4) pour un problème extérieur, comme le montrent les 
puissances de 7 qui y figurent. Pour le problème qui nous occupe 
(océan compris entre deux sphères), les séries devraient comprendre 
des termes des deux espèces, ainsi qu'il arrive dans le problème lapla- 
cien, sans rotation. 

Deux propriétés des fonctions ‘Tr les rendent utiles dans le pro- 


blème elliptique (w < - : 

1) Chaque fonction Trt permet de satisfaire à toutes les conditions 
internes et frontières, et donne ainsi une solution particulière. 

2) Ces solutions particulières sont orfhogonales, ce qui permet de 
trouver facilement les coefficients qu’il faut leur attribuer, pour obte- 
nir, par leur addition, la solution qui correspond à un potentiel per- 
turbateur Pe(x, )eïf quelconque. Ces deux propriétés font défaut 
aux fonctions Mi. 


ILE PRNNPRIÉIÉS US, ÉINILIIOS NPA MIISIIE 
1) L'ensesdie des deux termes 
2,74 LS ét 77) us 


D permet certainement pas de satidaise aux deux conditions fron- 
_ fésr=Rr=R — h, pes un simple choix des deux constuntes A, 
me en fonction de cée qi Ééfinit le terme correspondant de Pr, 
| parce que p see encore dans les deux équations (quelque simplifica- 
| fon qu'on fsxzne dans La forme des conditions frontière) en raison de 
Le différence de forme de MW, et de Mt_ys3, Onle reconnaîftra sans 
 Éfcdié, suisse des longueurs d'écriture, en essayant divers cas 
Dee sie par exemple Le plus simple de tous 

Pa? r=hk 

P=-0 1=R—h 


er sans se priéoccnges de La sipmification physique. P- doit étre déve- 
… boypéen TMS (sans M 4,1) 

_ 2) On s'assure facilement en exasminant La forme de M,, que l'on 
…— Na fpes sus le deni-mésidien entier de pôle à pôle 


Le. 


4 | 
Ÿ MG Mode = 0. 


| Pour écrise mue condition d'orthogonalité, 51 faudrait multiplier 
gas vu four d'imborians, où poids, je Eu 
des cite decmon, dei canine dus 
sales À, ». 

Kesistons plait sus La difiésence entre les M, et les M_,,,,, 
sa 

FPT, gs {3 — gp) Le 20) 
- TH = PER 2 — ps — pi? EF" 

On en déduit que les solutions (>) se composent de polynômes en y, 
‘avec le facteus 5 — 4% si hest impair, tout RL ue 
sales oréénrires. En ce 74) —%— es un polynôme eny,ouen 
nn come rh cu datée 
mais ce terne disparsit dans (:) avec le coefficient qui le précède 


L 


d Indiquons Les premières fonctions prises parmi celles qui sont 
1 PRE 


Su 
face à : 
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indépendantes de la longitude, et généralisent les polynômes de: 
LEGENDRE : 


MORT 

He 

mer [( er) Es] 
nm mn ne alt 


Les solutions (4), au contraire, qui sont : 


( 2n+i ; à ( 1 
May) = Ma X (sa) (ee HE) :) 


deviennent toutes infinies pour 4 — p., et sont à rejeter pour tout pro 
blème concernant la surface entière de la sphère. 

Voici où nous en sommes : 

Si la sphère est entièrement liquide, sans fond, chaque fonction 


1ka 104 
Aire M, u)e 


fournit la solution d’un problème particulier dû à un potentiel pertur-- 
bateur 


Pire M, x ( ue” 


si l'équation à la surface libre est P — Pre, 
Si l'équation à la surface est de la forme plus générale où apparais- 


L nus Lu AD x : » 
sent aussi les dérivées su” il faudra former les équations de récurrence 
Le 


entre les M et les . à l’aide de celles des T ; mais le défaut d’orthogo- 


nalité des M en rendra l’utilisation difficile, et fera probablement 
intervenir les parallèles critiques 14. 

Pour le fond sphérique, la solution paraît impuissante. 

On ne voit même pas quelle forme de fond équatorial (1? < 11?) 
ou polaire (u? > 1?) excluant les parallèles critiques, et indépendant 
de n, Rk, pourrait être traitée à l’aide de la solution en 7-(#+1), 

I1 faut donc chercher autre chose. 


5. — Pour écrire facilement les conditions aux frontières, toutes 
deux sphériques, il suffirait que ces deux sphères, R et R — x, fissent. 
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partie de l’une des familles de surfaces coordonnées, sans que toutes 
les surfaces de la même famille fussent sphériques, comme dans l’essai 
précédent. On pourrait donc chercher des solutions de l'équation hyper- 
bolique à l’aide de coordonnées curvilignes formées de trois familles de 
surfaces dont l’une comprendrait les deux sphères frontières. Mais, il 
se préserte, dès que les coordonnées ne sont pas orthogonales, des 
particularités qu'il vaut mieux mettre en évidence sur un exemple plus 
simple, celui du globe sphérique entièrement liquide. Je vais étudier 
une solution toute différente de la précédente, dans laquelle figure une 
famille d’ellipsoïdes, comprenant la sphère de surface libre. 
Reprenons les variables 


12 
HIME = 


qui transforment l'équation hyperbolique en équation laplacienne. 
L’équation de LAPLACE en x, y, £ admet parmi ses solutions des 

familles bien connues d’ellipsoïdes homofocaux de révolution autour 

de oz. Si l’on revient à z, on trouve qu’on peut prendre des coordon- 


nées #, v dans un plan méridien, définies par : 


( m—AV(T— #7) (1 — 1) 


z = Acur, 


() 


æ étant comme précédemment le rayon cylindrique. 
Les courbes # constant sont les coniques : 


2 


w° 2 de 
(2) A?(r — uw?) Arc — 


On obtient un cercle pour 
A?(r = U1°) = A?c°u;:° — RE 


ce qui donne 


I () 
Hu ROHA TEE =AVI EE 
C’est la sphère frontière cherchée. Si on cherche la famille des 
courbes v constant, la symétrie des formules (1) montre qu'on re- 
tombera sur l'équation (2). I1 est impossible de séparer, comme dans 
le cas del’équation, de LAPLACE les deux familles de lignes coordonnées 
en une famille d’hyperboles et une famille d’ellipses, deux courbes 
de la même famille n'ayant aucun point commun. 


MARCEÏI BRILLOUIN 


Si nous cherchons alors à satisfaire à notre équation par une expres- 
sion de la forme : 


nous obtenons : 


I _ EV | 
dv 1 — v° ( 


dv | RU ) + 


(| 4 
pipes, O1 — 46 | ) dv 
Nous savons d'avance que nous pouvons obtenir deux équations 
différentielles séparées pour U et V. Il suffit d'ajouter et retrancher 
n(n + 1)UV, # étant une constante que nous prendrons entière, et 
de diviser par UV. On obtient ainsi pour U : " 
d' n 4U RU ee 
ae —) %| a T6 + in + I)U — 6 JE 
c’est l'équation connue des fonctions sphériques tessérales : U — Tr(w). 


On obtiendrait de même la solution V — Tr“(v). 
En outre, sur la sphère u — uw, on a 


2.= ACHiv = R9; 
z 
R 


On aura donc une solution particulière de l'Equation hyperbolique 
en preñant 


et la variable v y devient égale au sinus de la latitude 


Pi — A he TK (u)T,A(v), 
qui donnera sur la sphère #, 
k — A KT & RO x +) 
P, ART, tu, )e TT, (& , 
où les fonctions Tnt (£) sont orthogonales, avec tout l'avantage qui en 


résulte pour l'identification d'une série S-P;* avec le potentiel trou- 
blant Pe donné. 


| ÉCHECS : 
L'avantage se conserve pour les dérivées +. ; mais avec ces 
Le € 
coordonnées ,v qui ne se coupent plus à angle droit, comment forme- 
rons-nous la dérivée normale, qui intervient nécessairement dans 
l’ensemble de la solution ? 
Etudions donc de près la définition d’un point au moyen des coor- 


données #, v. La première question à examiner est évidemment celle 
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de la réalité de #. Ce n’est en effet que pour # réel que nous pourr>1 
employer aisément les fonctions tessérales. 
L'équation des courbes # constant, résolue en #?, montre de suite 
que #? est réel pour : 
(em — 2 — cA) (cm — 2 + cA) (cs + 2 — cÂ) (co +2 +cA)>0 


ce sont là les premiers membres des équations de 4 caractéristiques 


CA É 

A PT 

à cÀ DR Ses 

\ 2 ÉS 

SE C ÉS 

ESA \ LL 1? LD KSS 
SOS , PP SL | LS 
Se [ CZ PS ASS 
SIA | ce CL | KS 
SSS | 4 Ç CD) LS 
ISSN ZX 


N 


(fig. 2) formant un losange qui a les axes pour diagonales ; #? est réel 
dans les portions ombrées du plan. Il est positif pour 


{ui — A?) — 2 0, 
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donc à l’intérieur d’une hyperbole ayant ses sommets au sommet 
du losange, ce qui exclut la portion quadrillée du plan. 

Il est clair que les coniques admettent pour enveloppes ces 4 carac- 
téristiques (1). La circonférence en particulier leur est tangente aux 
points 


! C 
SE ne) 
\ LUC 
c? 
(ne —- 3° 
TEFrC 


Nous retrouvons les points que nous avait donnés l’examen préalable 
des conditions de détermination de la solution. Ils correspondent aux 
«parallèles critiques » de la théorie des marées. 

; 4w? — 0? 

Si nous nous rappelons que c* — 7 -, nous voyons que pour, 
les marées semi-diurnes, c est très petit. Comme il représente le coefi- 
cient angulaire des caractéristiques, les latitudes critiques sont dans 
les régions polaires. Elles pourront n'être pas gênantes dans des pro- 
blèmes qui excluraient ces régions, ce quiest physiquement acceptable. 
Par contre, pour les marées diurnes, c est voisin de V3, les latitudes 
critiques sont d'environ 309, et les difficultés que nous allons exposer 
ne peuvent être tournées. 

Seuls nous intéressent les points intérieurs au losange. En ce cas 
toutes les coniques sont des ellipses. Celles qui correspondent à # > # 
(ou v > v,) ont leur grand axe suivant oz, les autres suivant ow. Un 
point intérieur au cercle est obtenu par l'intersection d’une ellipse 
d’une espèce avec une ellipse de l’autre, et correspond à w  w, < v 
(ou v < w, << u). On obtient un point extérieur, dans la région polaire 
en prenant # > w, et v > %,, dans la région équatoriale en prenant 
u Lu et 0 Lu. Mais au voisinage du point de contact avec le cercle, 
il est impossible de « sortir » du cercle. T,es coordonnées w, v ne se prêtent 
pas à la représentation d’un point extérieur le long d’un parallèle cri- 
tique. 

Les points intérieurs eux-mêmes donnent lieu à des difficultés. Au 
voisinage du point de contact les ellipses sont près de leur enveloppe. 
Ia profondeur de leur point d’intersection est petite du second ordre 


(1) Cette propriété correspond d’ailleurs à celles d’être homofocales pour les coniques fournies 
par l'équation de Laplace (tangentes à un quadrilatère isotrope). 


” 
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“par rapport à la variation de # (ou de v) à partir de #,. Pour obtenir 

par exemple une variation À suivant le rayon, il faut faire varier # et v 
SERRE h Jr 4? Heu 

de quantités équivalentes ren =, et — V A > : der 

Partout ailleurs une variation d7 du rayon vecteur est linéaire en du, 


ù HAE. Mes dP : ; ; 
d'u : la dérivée + est linéaire en Mais à la latitude de contact 


(v — %), une variation dr est proportionnelle à (4 — w,)?, (v — #,)?. 


À cette latitude (Note IV), la dérivée Æ devient linéaire en = et 


\g 


Fis: 3. 


2° ST z PAPE : . 2e 
>: De là résulte la nécessité d’une discussion spéciale de toute la 


solution et de la condition frontière, au voisinage de cette latitude. 

I,a même circonstance complique l'écriture de la condition au fond, 
quand on essaie d’utiliser cette solution pour l'étude (approximative) 
d’un océan de faible profondeur. Loin des latitudes critiques, une 
faible profondeur est proportionnelle à (4 — z,) ; aux latitudes critiques 
elle est proportionnelle à (# — ",)? ; cela exige des précautions dans 
l’approximation. 
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6. —— Pour terminer, nous allons montrer que nous sommes loir 
d’avoir épuisé toutes les formes possibles de développements en fonc- 
tions simples. 

Partons pour simplifier de l’équation de LAPLACE 

NEO; 
et supposons que nous écrivions la solution P sous forme d’un produit 
de deux fonctions. Pour voir à quelles équations on peut assujettir 
celles-ci, posons : 
» à à” 2° o® ? 
P— 6c _ d'où ue (her 
oX 2% \dX , 


L'équation en ® est : 
A ct A, A 0, 
en posant 
do \? dP \? o® \? 
te EE 
La division de P en produit se traduit par le partage de ® en somme, 
soit : 
Lo = b, —- Ds. 
On a alors : 


do; db: à d®P; d®;, re CLR ee à 


d D, + Ad, = — 
(3) Adi + AŸ, AP; + Aiès Rs dX dy 2 d2 22, 


Toutes les solutions qu’on a l’habitude d'examiner sont telles que 
le second membre soit nul, c’est-à-dire telles que les surfaces ®, cons- 
tant, ®, constant soient orthogonales D, et ®, sont alors des solu- 
tions de : 

| Ab, + A,d, — A 
| Ab, + 4%, — — A, 


(A étant une constante arbitraire), assujetties à la condition simple 
précédente. 

Mais on pourrait aussi rechercher des solutions telles que le second 
membre de (3) soit la somme d’une expression fonctionnelle dépendant 
de ®, et d’une expression dépendant de ®,. 


( 2 
AF,(P, e is 0) 
OX TOUR 
] + AS op, 0, E db; ac 
2 »Y b] A ? 2 , , 
oX” à ox 


= 278 + + 
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Ces expressions pourraient être des fonctions de ®, (ou ®,), de ses 
diverses dérivées, des variables indépendantes, comprendre au besoin 
des intégrales. 

Par exemple si l’on prenait : 


| ALL 2P; oP, 2P; dP4 2e A d°®, he do, c o°P, 
(4) \ ax 4 dy 2y 22 02 ) DECO E RETRAITE 
do? do? o°®, 

A ET RRQ Er 

( D ETS ? 27? CG 22? ? 


A4, B1, Ci, Av, Bo, C étant des constantes, ®, et D, satisferaient aux 
équations : 


ee 


(x + A)” Pal Me 


—- A,®, — À 


sine 


raie + (I D. Han) + Ad — — À, 


qui ne sont guère plus a que les HAE La dificulté 
consisterait, bien entendu, à choisir des solutions qui vérifient la condi- 
hion (4). Il sembie probable que pour certaines formes de condition 
remplaçant la condition d’orthogonalité on aurait encore des résultats 
simples. Leur étude systématique présenterait le plus grand intérêt. 

Les mêmes idées sont applicables à notre équation hyperbolique. 
Voyons seulement par quoi est remplacée la condition d’orthogonalité. 


Si on change zen £ dans la condition d’orthogonalité, elle donne : 


oP, oP; Po o®P; c 9P, Ps 
dX 0% dy ày D201021 1. 


Supposons que l’on cherche des solutions qui soient de révolution 
autour de 02 ; la condition se réduit à : 
ob, 2P Lee a, do; A 
D 05 dZ 02 
Elle exprime qu’en chaque point les courbes %, constant, ®, constant 
et les caractéristiques forment un faisceau harmonique (). En parti: 
culier si l’une des courbes devient tangente à une caractéristique, 
l’autre lui est tangente au même point. On retrouve, généralisée, 
une des propritéés de notre dernier développement en série. 
Quant à la condition d’orthogonalité vraie, satisfaite par les termes 


(x) Il en serait évidemment de même dans toute transformation homographique. 
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.de notre premier développement, c’est ici un exemple de condition 
plus compliquée imposée aux fonctions D, et P,. On vient de voir que 
c’est un cas particulier parmi l’infinité des conditions possibles. 


7.-—Mon but, dans ces conférences a été d'attirer l'attention des jeunes 
mathématiciens sur quelques difficultés, intéressantes en elles-mêmes, 
et plus encore pour les applications hydrodynamiques, qui se présentent 
dans l'étude de l'équation hyperbolique la plus simple, dès que la 
frontière n’est plus celle que les mathématiciens ont pris l'habitude 
d'étudier. J'espère que ces difficultés exciteront la curiosité de quelques- 
uns d’entre eux. 


Note I.— Dès la fin du xvrrre siècle, les géomètres qui se sont occupés 
de la théorie du son ont rencontré des équations Avperboliques, et se sont oc- 
cupés de leur « intégration générale » à l’aide de fonctions arbitraires, ou 
de la détermination de solutions particulières convenant à des conditions 
aux limites données (problèmes pour lesquels l'intégrale générale fut re- 
connue être rarement efficace). On pourra parcourir avec intérêt le vieux 
traité de calcul intégral de Léonard EULER, les Œuvres de POISSON, certains 
mémoires de Du Bois-REYMOND (1864 ; 1883) et enfin de RIEMANNN (1860). 
Leurs résultats généraux se trouvent résumés, par exemple, dans le Traité 
d'Analyse» d'E. Goursar, t. III et les « Leçons sur quelques types simples 
d'équations aux dérivées partielles » de E. PICARD (Cahiers scientifiques, 
1927). Le cas de plusieurs fonctions variables et plusieurs équations aux 
dérivées partielles a Eté étudié par VOLTERRA (Membranes), | Acta Math. 
1894]. 


Note II. — Considérons donc le mouvement relatif d'un fluide par rapport 
à des axes oxyz tournant autour de oz avec une vitesse uniforme w. En 
variables d'EULER, les équations du mouvement sont du type : 


ou ou ou ou Mo 


— +u— +v- D — = X — —<. 
ot dX SE QU A d2 dx 


X, V, Z représente la force totale, comprenant la force de CORIOLIS. — 2wv, 
2wu, 0, et les forces attractives dépendant d'un potentiel V. Nous pose- 
‘TONS : 

P=V—ph, 
-et nous supposerons les vitesses assez faibles pour que les termes du second 
ordre soient négligeables. Les équations se réduisent alors à 


ou oP 
ss 2 —_— E 
\ ot + ox 
LS ov oP 
‘(x) { : — QU —= > 
| ot dV 
ow oP 
ot _ d7 
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Si on élimine #, v w entre ces éq 1ations et l’é juation de incompressibilité,. 


dU dy , ow 


—— 


PLAT ON VCSr 0 
on obtient l'équation en P : 
o°P 


d° \ 2 
ft (AP) + 4° 32 


= ©. 


Supposons alors que P, #, v, w soient des fonctions périodiques du temps, 
que nous écrirons (en changeant la signfication des lettres), Pet, weidt, velbi, 


; £ = 14 : 
we”. Nous obtenons enfin l'équation : 


®P . 


dP .4w? — 6? 52P 
dx? d "y 


oy? 42 22° 


(A) 


donnée dans le texte (!). 

Précisons les conditions du problème. Le liquide est supposé recouvrir- 
un globe sphérique. La vitesse de rotation étant faible, nous supposons que 
la surface du liquide tournant et gravitant ne s’écarte pas sensiblement d’une 
sphère concentrique au fond. Des forces perturbatrices données dépendant 


du potentiel P“ei5: sont alors introduites, imposant (?) la période Fe. Nous 


négligerons l'attraction exercée sur le liquide par le bourrelet (positif à cer- 
tains endroits, négatif à d’autres) qui se forme alors. Les forces sont donc 
toutes connues. Il reste à écrire les conditions aux frontières, savoir l’uni- 
formité de la pression à la surface libre, et la condition au fond. Il est com- 
mode pour cela d'introduire les coordonnées polaires : 7 (rayon), y sinus 
de la latitude), + (longitude). Klles permettent d'écrire la composante de la. 
vitesse suivant le rayon sous la forme : 

20 


U, — 


2 


92(4w* -— 1?) 


qui s'obtient aisément à partir des équations (x). 
La condition à la surface libre donne alors (compte tenu de la condition 
correspondante avant l'établissement du potentiel P°) : 


(D — P) + gUr = 0, 


toutes les quantités étant prises à la surface libre de repos : 7 — KR. Dans- 
cette formule, il faut choisir convenablement l'intensité de la pesanteur g, 
en fonction de la latitude. 

La condition au fond s'écrit simplement : 


De 10; 


h étant la profondeur du liquide. 


(x) C£. Poincaré, Acta Mathematica, t. VIL (1885), p. 356. 
(1) Au bout d’un temps suffisamment long, en vertu des résistances passives. 
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Enfin, si l’on veut avoir le soulèvement « en chaque point de la surface, ce 
qui est le but ordinaire des problèmes de marées, il suffira, P étant supposé 
-connu, de remarquer que : 


Ainsi tout revient à déterminer une solution P de l'équation (2), régu- 
lière entre deux sphères concentriques, et satisfaisant sur chacune d'elles à 
“une condition linéaire par rapport à P et àses dérivées partielles. C’est ur 
problème classique dans le cas de l'équation de Laplace à laquelle se réduit l'équa- 
tion (2) pour w — 0. Mais l'équation (2) est une équation hyperbolique, 
au moins dans tous les cas suggérés par les marées terrestres. Pour les marées 
semi-diurnes luni-solaires, 0 est en effet très voisin de w, mais plus petit. 
Pour les marées diurnes 1 est voisin de w. A fortiori, pour les marées semi- 
mensuelles, mensuelles, etc., le coefficient de, est toujours positif. 

Dans ce cas, rien de classique ne permet de former la solution cherchée. 
Dans les problèmes étudiés de propagation du son, la variable f remplace 
la variable z, et rien d’analogue à la distance, aucune expression de la forme 
x? + y? + cl, ne saurait intervenir. Danslestravaux de POINCARÉ sur la 
même équation (1), la nature de l'équation n'intervient pas, elle est traitée 
partout comme elliptique. Enfin la forme sphérique des frontières n’ap- 
porte pas de simplification importante, comme dans le cas de l'équation de 
LAPLACE, laquelle possède, peut-on dire, la symétrie sphérique. 


Note III. — Posant 


PA F 
- = b, =. 4?) + y? Li 


LÀ 
É SANS ON , o2M : RE à 
(0? + 4w?pu?— 4w?) (x — 1?) us + [w?(12 — 8n) (1 — pu?) — 20% ]u Er | 
np 6! 
+ [1308 — qui] + nm — 2) (08 — gruta) — ŸE [n  À 


l 


Remarque : Une fonction M* ne convient qu'à une seule valeur de », 
puisque # apparaît dans deux coefficients — et linéairement dans l’un d'eux. 
L'équation correspondante du problème elliptique (fonctions tessé- 
rales T,*) ne contient n que pour leproduit # (n + 1), qui ne change pas 
si on remplace # par — (n + 1), de sorte que l'on a T,f = 'T#_,,,, 


Note IV. —— Des formules 


o = AV — wi G — wi), 


Ly1 
| 


— Acuv 


«(1) Loc. cit. 


sn HER 


QUELQUES PROPRIÉTÉS D'UNE ÉQUATION HYPERBOLIQUE 


on tire 


du — cUSdS + Av(r — u')dz 
cA?(v? — u?) 
_cvwdw + Au(r — v?)dz 
CA (0? — y‘) 
= Ar — un — 0° + (1 + c'hu?v?] 


dv — 


ee 
Æt, à = constant 


AN UT MR Te 
3 (I ee (I ve 


d (02 — w?)A VI — ut — v? + (1 + PER 


Le dénominateur v?-u? crée évidemment une difficulté pour v = w, c'est-à- 
dire au point de contact avec le losange des caractéristiques. 


En particulier, pour la solution T,A{u)-T,*(v), | et | deviennent égales, 


, , Prier 
pour # = v ; de sorte que le numérateur s’annule aussi ; É prend un as- 


pect indérerminé. Mais si, pour la même valeur de #, par exemple #,, on 
développe en (v — w) le numérateur, on obtient sans difficulté 

aP 

dv? 


— 2 >P Re 
lot 
W;v—=u A-QuV1 — 20? + (x + c'jut 


(Deux conférences données à l’Institut H. Poincaré le 28 et le 
31 Mai 1930 et rédigées par M. Coulonel). 
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Sur quelques problèmes non résolus 
de la Physique Mathématique classique 
Propagation de la fusion 


PAR 


MarcEz BRILLOUIN 


I 


1. — La Physique mathématique classique, celle du x1X® siècle, pose 
des problèmes de deux espèces très différentes : 

19 Ceux dans lesquels les domaines dont on étudie le changement 
de propriétés sont limités par des frontières géométriques données ; 

29 Ceux dans lesquels ce sont ces frontières, soit fixes, soit mobiles, 
qui sont les véritables inconnues de la question. 

Pour les problèmes de la première sorte, on est en possession de 
méthodes mathématiques variées qui conduisent généralement à la 
solution requise. Mais pour les seconds, c’est presque aussitôt après 
la mise en équations qu’on est arrêté, et dépourvu de toute méthode, 
ou même d'exemples à imiter. 

I1 suffit de rappeler combien peu d'exemples on a réussi à traiter de 
la forme des gouttes liquides ; quels détours KIRCHHOFF a dû em- 
ployer pour l'étude des courbes élastiques. 


2. — HÉLMHOLIZ ayant indiqué la condition d'existence de sur- 
faces de discontinuité permanentes dans les liquides sans viscosité, a 
formé quelques exemples : KIRCHHOFF a donné d’abord un peu plus de 
généralité à la méthode mathématique, puis l’a bien présentée, mais, 
pour les mouvements plans seulement, grâce à des propriétés par- 
ticulières des fonctions imaginaires. La méthode réussit parce qu’en 
permutant les rôles des variables géométrique et dynamique, on est 
ramené en fait à traiter un problème à frontières fixes données — 
grâce à une représentation conforme. L’intéressant ce [serait l'es- 
pace à 3 dimensions, et les surfaces de discontinuité qui naissent des, 
bords d'obstacles limités en tous sens. Ici la permutation des variables 


géométriques et dynamiques ne sert plus à rien. 
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Sans discuter ici aucune des questions d'ordre physique que sou- 
lève le problème hydrodynamique initial, je me borne à constater 
l'impuissance en présence du problème mathématique précis que voici : 

Etant donné un obstacle, à bords aigus, plongé dans un courant 
liquide uniforme et constant en amont au loin, trouver la surface 
de discontinuité permanente émanée des bords aigus, le long de la- 
quelle la vitesse tangentielle est sur la face amont égale à la vitesse 
du courant, et nulle sur la face aval, où tout le liquide est en repos. 

Arrêté dès ce premier pas, le géomètre reste coi en présence du pro- 
blème d'état variable correspondant et de celui des pales d’hélice, 
même schématisé au delà de toute vraisemblance. 

Au technicien de se débrouiller comme il pourra. C’est une situation 


peu glorieuse. 


8. — Je prends une autre question posée depuis FOURIER, — un 
siècle — la propagation de la solidification du globe terrestre, supposé 
en fusion au temps origine, et elle aussi, schématisée au delà de toute 
vraisemblance, au point de faire perdre à la solution toute utilité géo- 
désique, sauf comme tout premier travail d'approche en vue de pro- 
blèmes un peu moins éloignés des conditions probables. 

C’est ce premier problème que je vais poser, et dont je vais mettre 
en évidence les difficultés : 

Propagation de la fusion de la glace. — Propagation de la solidifica- 
tion d’une croûte sphérique. 

Le problème plan le plus simple a été traité il y a un siècle par 
LAMÉ et CLAPEYRON, dans un court mémoire, lu à l’Académie des 
Sciences le 10 mai 1830 (!). 

Dès que la température initiale n'est pas uniforme ; ou la tem- 
pérature superficielle pas constante ; ou la profondeur totale limitée, 
les difficultés semblent insurmontables. 

Et pour les problèmes sphériques, c’est bien pire. 

(1) Publié aux Annales de Chimie et de Physique t. XLVNII, 1831, p. 250-256. 

Malgré son titre « Mémoire sur la solidification par refroidissement d’un globe liquide », c’est 
seulement en paroles qu'il s’agit d’un globe. Le problème réellement traité dans la note finale 
de deux pages, est celui de la formation d’une croûte solide plane par refroidissement de la 


surface supérieure à une température constante, le liquide indéfini au-dessous restant à la 
température constante et uniforme de solidification. 


s = BVi. 


L'épaisseur croît proportionnellement à la racine carrée du temps ; mais les auteurs ne déter- 
minent pas le coefficient 8. 
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II 


4. — Surfaces 1sothermes planes. — Commençons par un problème plan 
dont Franz NEUMANN (1798-1895) a fourni la solution dans ses Leçons 
inédites de Künigsberg au début de sa sixième année de cours (vers 
1835-40 ?) d’après H. WEBER (1). 

La même solution a été donnée par STEFAN, d’une manière indépen- 
dante, en 1889, à l’Académie des Sciences de Vienne, et imprimée dans 
le volume XCVIII des Siézungsberichte de cette académie (Abth, Ile) 
P. 473-484 paru en 18900. 

Soient : 

c; la chaleur spécifique de la glace ; 

Pa Sa densité ; 

À, sa conductibilité calorique ; 

Cas Pa R2, les mêmes grandeurs pour l’eau ; 

L, la chaleur de fusion de la glace ; 

©; la température de fusion ; 

z la profondeur depuis la surface libre de référence (z = O) ; 

@, (z, t) la température sur le plan z dans la glace, à l’époque ? ; 

O, (z, t) la température dans l’eau. 

Pour réduire le problème mathématique à ses difficultés essentielles, 
on le simplifie en laissant de côté des propriétés physiques pourtant 
importantes : 

Simpli fications. 

On suppose l’eau liquide immobile maintenue à la température 
®, > (}; constante, pour z très grand. On néglige les variations de 
chaleur spécifique, et, ce qui est plus grave, de conductibilité, avec 
la température. On néglige également la dilatation, c’est-à-dire les 
variations de densité. 

En résumé on traite ci, R1, 1, C, 02, kx, Comme des constantes. 

On néglige aussi le changement de densité par fusion. Toute cause de 
changement de pression ayant été négligée, on peut supposer que la 
température de fusion @/ est une constante, ainsi que la chaleur de 
fusion I. 


(x) Die partiellen Dijferential-Gleichungen der Mathematischen Physik, nach Riemann's Vor- 
lesungen, Braunwschweig, 1901, 4° édition, t. II, pp. 118-122. 
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5. — Le problème ayant été ainsi schématisé, les équations sont 
or 
d2? of 


\ 2°@> __ 28e 
SAUT 0e 


\ dans la glace (À) 


| dans l’eau ( ke 
| CeQ 


\ 


à la surface de jonction z = €(é) 


RE 0; — 8, — 8; 
avec 
où {, 20, 
! — sé Rues ba. 
IX Ep = (A 22 Re St 


On écrit ici s (sans distinguer entre la densité liquide ou solide l), 


ns 
dans le 1° membre, où (0%) est l'épaisseur d’eau congelée par se- 
conde. 

Le second membre représente, l'excès de la chaleur enlevée par 
conductibilité du côté « glace » sur la chaleur apportée du côté «eau ». 

Même dans ces conditions physiquement simplifiées à l'excès, on 
pourrait être curieux de savoir comment l'épaisseur de glace varie 
avec le temps, pour une distribution initiale, @, (2), arbitraire des 
températures dans toute la profondeur (2 = o) avec une loi de varia- 
tion Ok (é) arbitraire de la température superficielle. 

Ce problème n'a pas été traité. 


6. — Le problème traité, parce que des intégrales de l'équation de 
la chaleur déjà données par FOURIER et LAPLACE vers 1820 (1}, s'y 
sont trouvées adaptées, est le suivant. 

Etat inihral : pas de glace, © — 0 ; température initiale de l’eau 
uniforme Oj(z) — @,, quelconque supérieure à la température de 
fusion. 

Au temps é — 0, on établit brusquement la température @ < €} 
à la surface libre z2 — 0, et on la maintient consfante en absorbant 
toute la chaleur que la solidification de la glace apporte par conduc- 
tibilité à la surface libre. 


Aïnsi, la solution s'applique au cas où, au-dessus d’un lac profond, 


1) Voir note finale. 
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dont toute l’eau est à une même température quelconque (+ 5°, 
+ 109, + 159, par exemple), arrive brusquement de l’air au-dessous de 
zéro dont la température se maintient exactement fixe (— 10°, par 
exemple). 

Il s’est formé instantanément une mince couche de glace, dont la 
face libre est à O4(— ro°) et dont la face en contact avec l’eau est 
remontée aussitôt à la température de fusion @. 


7. — La solution est 
2 19 

ér1( V “) 

ërf(> V je) 


| (EVE) ef 


| O1(2; 1) — 8; — [8 — 8,]| 1 — 


III 


8:(2, Î) ra 0; Eu ke 0;] 


en posant, suivant une notation qui rappelle qu’il s’agit de la fonc- 
tion des erreurs. 


u { &rf(o) = 0 
(IIT') 8rftu) = ? JL CR 0 DR 
VT/o | &rfloo )—7x 
La surface de contact eau glace, est au point 
a doit être tiré de l’équation 
| + a? cp 
2eù VITE che 22 eg Mar k 
FM Te cp 
IV a? C0 “4 e 2 
Rap ec 6, + | 
T * à C20 
I — &rf é \ ) 


Pour que la solidification progresse vers les z positifs, il faut que 
les températures extrêmes @%:, ©, soient telles que l'équation donne 
une valeur de x positive. 
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La discussion de l'équation en x n'a été entreprise ni par Neumann, 
ni par Stefan, ni par H. Weber. 

Pour une matière donnée il y figure 5 grandeurs caractéristiques de 
cette matière : 


cp CP 
1 ; ko n ka, ka, RLe. 


I1 faudrait dresser le tableau numérique des x, en fonction des 
deux données qui restent arbitraires : 
8,—8,<<0, 8 —8,; > 0. 
On peut adopter pour l’eau, comme ordre de grandeur (CGS) 


9,90) ne 
0,004? 0,00125 ? 


0,004, 0,00125, 80. 


\ 
L'équation en ( ; } est alors 


0,045 (8; Et 6,) £ 


CTES :) 


1 —6(28 ©) 


En se donnant un certain nombre de valeurs de x, on obtiendra un 
nombre égal de relations linéaires entre @ — © et @., — 0x. 

On tracera les droites correspondantes, avec leur cote %. Si le ré- 
seau ainsi tracé est assez serré, on déterminera à vue, d’après la posi- 
tion du point donné @; —@;, @x — A} pour chaque problème parti- 
culier — la cote de la droite du faisceau sur laquelle il se trouverait. 
D'où le coefficient de V/# dans la propagation de la fusion. 

Répartition des droites : Calculs et tracés de M. Couroms. 


— 0,035 (80 Fr 8;) 


RE — 800% 
45 —7— (8 — 0) — 35, (dx — 8) = 142000 À 
6rf (xx : T — 6rf (28 :) : 

ou 
AX — BV = C 
I — «x petit. 
SI = 27 
1) SR ESS 
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Si on se contente d’abord du premier ordre d’approximation, pour 
voir l’allure du faisceau, on obtient les droites 


3:63 ÿ 
(2) 


qui enveloppent la parabole V — — 13/— X (fig. 1). 

Exemple : Pour 4 — 0,01 on obtient 726 X — 35 V — 710 au lieu 
de 726 X — 40,6 Y — 710 (voir tableau I). 

Une bonne approximation est obtenue par la formule 


— 35 Y — 142000 CG} 


303 X — 35 (x m0 :) Y — 142000 = 


(2) 


où 2 erreurs de signe contraire sont commises sur le coefficient de Y. 
Exemple : Elle donne pour # — 0,01 : 726 X — 40,5 V = 710. 


YŸ |la=0 (À pett ïl 


X x 
a er 
Fig. «. one 
II. — à grand 
re I 
&rf(x) = 1 — 2. \ 22 SE ) 


On obtient de même au premier ordre d’approximation 
a? a? 
— 125 — 7 
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ou sensiblement 
2 
6 [e 2 


Ve — ge" es 


La droite s'éloigne rapidement à l'infini en devenant parallèle à ox 


(fig. 2). ° 
Exemple : « — 0,4 donne (par le calcul complet) : 


0,328 X — 224 V = 56 800. 


On voit avec quelle rapidité la droite s'éloigne. Ce cas n'offrant 
aucun intérêt pratique, nous ne chercherons pas de formule d’approxi- 
mation. 

Tableau applicable à la glace : 

On peut restreindre à 0 € X << 20°, 0 < V < 800. D'où : 


Points ayant servi à la construction des graphiques. 


Y=0 NE Cotes en dixièmes 
dou X=1,0 |[X—5,5 de degré. 
V = 0 V — 80 
TT KReuo Rent 
sa NET NX 120 
OS] x= 00 | v = 489 
en NAT 0 NN —120 
Fr X = 162 | V = ri, 


C’est à l’aide de ces nombres qu'a été construit le graphique. Dans 
leur partie utile, les droites sont tracées en trait plein : leur prolon- 
gement, en traits interrompus, a servi à tracer l'enveloppe. Celle-ci 
est assez bien définie pour permettre de tracer sans erreur notable, 
la droite qui passe par un point quelconque @; — OO, 2, — ®}, 
donné et de la coter en 7, par interpolation à vue (par exemple sur 


s 


SUR QUELQUES PROBLÈMES DE LA PHYSIQUE MATHÉMATIQUE 


Pig. 3. — Eau et glace : 6, — @; en ordonnée ; 0; — 6, en abscisse.- 
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l'arc de cercle tracé du point de contact avec l'enveloppe comme 
centre). Partant de cette valeur de %, l'équation fournira par approxi- 
mations très rapides la valeur exacte avec 3 ou 4 chiffres si on le juge 
utile. 


8. — STEFAN a résolu le problème un peu différent dans lequel une 
colonne de glace infinie du côté négatif, à température uniforme @ 
est subitement mise en contact avec une colonne d’eau infinie, à 
température uniforme ®. La face de contact prend subitement la 


température de fusion @;, la face extrême, z — — , de la glace est 
maintenue à (y, et la température ‘le l’eau au loin (2 = + ) 
Tee oO 


Les températures sont partout données par 


ENCORE TEEN C 
IE k) érf(2 Es) 


a /Ce 2 À ce 
| érf( ) —r(5 .) 
8 — 6; — (8, — 8,) - - D PER 


avec la propagation 
Vo) Éeavt 


et la condition thermique 


RE 
RE à (CP à AA _ 2, 4h 
/ Se / 
| I +8 Ve) t VTT 
2V k; 
(VI) : 
, = __ @°cop 
+R air = Li VE Up 4ke 
PEN 
\ k; / 
9. — Arrivons aux problèmes plans, avec distribution initiale 


“arbitraire des températures. 
Commençons par le problème général plan. 


" 
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Une couche de glace d’épaisseur Z, est surmontée d’une couche d’eau 
d'épaisseur 2. 

Dans l’état initial, la température de la glace est ®! (2), et celle de 
l’eau ®7 (2). A la surface plane de contact, la température commune 
est la température de fusion @j. 

On fait varier d’une manière connue en fonction du temps les tem- 
pératures @E (— Z,, #), On (+ L, #), des plans extrêmes z = — L,, 
Di Le | 

On cherche comment progresse la surface de fusion, et comment 
varie la distribution des températures. 

Voici une manière de poser le problème mathématique qui conduit à 
“un système de cinq équations intégrales : 

Posons 


tki\ er, C1£ NCA TT 
(VIN) N(z—, ea) = Vie 
L//équation de conductibilité a, comme on sait, (LAPLACE, FOURIER) 
pour intégrale générale 


(VIII) OZ D == f CN (z — Y, des ) 
CPE ee) C1P 

où ®! (z) représente la distribution initiale des températures @, (z, 0), 

tout le long de la colonne — © < 2 < 0. 

Pour une barre de glace (I) limitée, — /, < z < 0, la fonction ®! est 
donnée, et connue, seulement dans l'intervalle — Z, -> 0 ; en dehors de 
cet intervalle elle est inconnue. Ce sont les fonctions 
| ECUECEE STr 
(IX) ) : : 
L p— (y OS ZI co 
que je vais prendre comme inconnues ainsi que les fonctions ana- 
logues relatives à la colonne d’eau (II), qui s'étend de o à + Z. 


(eo) 


/ PE Ro: 
(X) On(z, é) = ICONE =} LE 


avec 1 donné le long de la colonne o < z < /, et inconnu en dehors 


RME UOTE 


SH NIO 


8 


(XI 
2 RER ee) LD E<ESto. 
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Pour déterminer ces quatre fonctions W ainsi que la loi de propagation 
(XII) z = {(t) 
de la surface de fusion, on aura les cng conditions ‘aux extrémités 
libres fixes, et à la surface de jonction mobile : 


Extrémités fixes : températures données en fonction du temps, de- 
puisi—= 0 
Fr 00 
QI(— /,, é) — | _ PON(— pe Nr 


NE 
(XII) pa 


On}, {) — | UN (4 — }, ms dr. 
— co 


Le 


Surface de jonction ; température de fusion Oj, constante : 


F t 
paf MN — x, a) dy 
(XIV) FE 


@, = |" enr Von the NL 


Equation ORRESR RE 


VE th: 
life fo DN'(E) — x, F2) dr 
(XV) 
| — ke : BICUN (LE ENT te) dy 
e/ — © Cap 
en posant 
LOI e tk 
N'— 2 2 - x N(: etes 


Les deux premières équations sie aux extrémités fixes, sont de 
types connu. Mais l'introduction de € (£) à la place de z dans les noyaux, 
N;, No, N;', N2', dans les trois équations de la surface de jonction, four- 
nit, pour la détermination des 4 fonctions WT, WI, Wr, Wn,et de Ë (é) 
un type d’équation intégrale, non étudié jusqu'ici, du moins à ma 
connaissance. 


Le problème est posé ; il serait bien désirable qu’on sût le résoudre. 


10. — Fssayons d'appliquer la méthode précédente au Problème de 
STEFAN. 


ETAT INITIAL 
Glace — 0 T7 < 0, O(z, 0) — 6, constante 


Ÿ 
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Alors : 
b, =:0,: Y, est inconnue, w_—0o 
Eau 0 < Zoo O(z, 0) — 8, constante 
Pi; Y_ est inconnue, Pate 0: 


L'extension infinie des deux colonnes supprime deux fonctions in- 


connues. 
EXTRÉMITÉS 


Z—=—@, 8, — 68 
: ! constantes 


2—= + ©, 05 =——= 6; 


À une époque { quelconque, la température de la glace sera donnée 


par 

2 je 
I [ AN 
= 6, 


8,(z, 1) == — 
Vz  — œ 
2 


et celle de l’eau par 


ee) 
I — 

8:(z, {) — —| 68. ( __€ TMdaB + k 
ia De RU: Eee 


Les conditions aux extrémités sont satisfaites d’elles-mêmes, ce qui 
supprime deux des équations de condition. Restent les deux conditions 


Le) EE 
22 Rit\ —62 
QT + _w;(z+28 SE Î ag 
JET MER Cap 


» 


thermométrique et calorimétrique à la surface mobile de contact : 


PR 2 a 2% 
LE 2 \ 1Ë 22 L 
vi Vr e, |. # He AtE \ 
2 kit 
5), /p 
. L. _&œ.. Ty Vite 
lo mr PT Lee 
x} 6. 59, /cp /— 00 : 
avec 
__E?ae co 
k| 8,—%,(0) ,— ant F 
nt Je, /& 
26 vr / a Si 2 V kit 
Les — 
{bee ,—tÙ V# 
ke 8, = w_{(o) —= ARot | 2 Rot 
— 7 = {/cpe + 
vz 2V ot U — 


Le lt) 
L 4 (et6)+28 L ag 
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On resterait encore bien embarrassé pour tirer € méthodiquement 
de ces équations, si leur forme ne suggérait d'essayer € — ay. et aussi 
— à tout hasard — de prendre Wï constante, et W1 constante, à 
choisir convenablement, s’il n’apparaît pas d’impossibilité. 

I1 faut alors 


I aV Cp ay C10 
re 2er — : r1( CU )| or va]: Le A | 


= eff) + 2e fe] 


Puis 
[ __ ce 
ke: FNaPre. — Ye n — 0 
ca VAE 
25 __ 63p 
fe Ra Ve: (e, = we # +- 0 
Vr Vhe 


en supprimant le facteur commun ï ; 
C’est précisément l'équation VI de STEFAN pour la détermination 
de x ($ 8). 


11. — Le problème de NEUMANN se présente un peu autrement 
parce que la colonne de matière est limitée au plan z — o, et s'étend 
seulement du côté positif. 

Pour une telle colonne de matière, si elle restait homogène, il fau- 
drait connaître la distribution initiale @ (z, 0) des températures sur la 
colonne 0 < 2 © et la loi de variation @ (0, é) des températures à 
l'extrémité en fonction du temps 0 < {<< %. 

On aurait alors, comme on sait 


BE SULER _s+y'æ 
re : "REA ykt Se 
GEAR NC CIE DRE ee 
(XVI) Vo 2VrTkRt 
+1 5 2co d 
+ [ e(o,t)e RG) AVR de er 
\ ./0 2VTR (4 ts 


Dans le problème de NEUMANXN je regarderai On (0, {) comme une 
fonction inconnue, pour l’eau. Pour la glace, au contraire, qui va naître’ 
à partir de l’origine, c’est la distribution initiale @: (2, 0) qui est in- 
connue, et la température à l’origine @' (0, {) qui est donnée. 
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Les conditions thermométrique et calorimétrique à la surface mo-- 
bile de jonction vont fournir les trois équations intégrales nécessaires 
pour déterminer ces deux fonctions @r (4, 0), OM (0, f) et la loi de pro- 
gression z — E(#) de la surface de solidification : 


Je poserai : 
@1(0,{)— Bt),  8I(z, 0) = w,(2) glace- 
O0, £) — Wu(f), OU(z, 0) — p(z) eau: 
Les ® sont données ; les Ÿ inconnues. 
1A/ cie co | 60 ae — 61 + vf 
OP Rte W1/-+)| e 4kit = 2 4kRit la 
HE Fee se GI Ÿ 
en ï4 LT dr 
“ a t C10 SR LT) LA 
+ 0,/2e [ Pre  4A1(é d'a 
XVII 
) 1, Je co | __{(6(9 — y)?cp el 
DERe Ne Ê ot ARE EC Lee 
A / © 
Et)” coo ä 
LE, / cp Tee AR À 
dl M er 
Puis 
ï Je 3/ So | (60e )?c1p 
Ne. ( ira [ VANLGE)—ne 4 
«/0 
| 
| eme Coco |s 
» JO 
CUT er 
Hi Lo) [ CE TON eme, 
A £() C0 
\ Ricio £ MED 2 £ - 
ot 2VT EP Ao) F) le 
/ D re 
(XVIII) Le — ee as | (K()— y)2c00 
ia nee 07 [ DCE (EE — ve re 
4VT Rat e/ oO 
E (£(4) + y}2c00 
— (bye ta 
E(#)*co0 
I 3/9 —2 s À TMEATES _& 
ER A pen 
= . C()”cop E 
2VT o ( =) É 
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Le cas général paraît encore moins abordable que dans le pro- 
“blème de STEFAN. 


12. — Examinons donc le cas où les fonctions données sont des cons- 
tantes. Tous calculs faits, il vient 


8, = dl, + (WI — 7 V #) 


et 
__WW°ap = _ {Dee 
Lo — VÉCOUX = Ptje 4kit Le a (wu te pu)e 4kot 
ro 
Avec 
&(£) = avé 


.ces trois équations deviennent 


or 1/%) + of: —wi(2 V5) | 
2— smér( © œ) + + val — erf(2 2 /Se)] 


-et, en supprimant le facteur Le) 


vi 
aci0 a°C9 


Loz PU gear SI __ RATES LT WT __ pu = 4ke 
ee = (WI — ie NE Vire (4 æue 
-.ce qui reproduit la solution même de NEUMANN ($ 7), (IV). 


13. — On peut se demander s’il n’y aurait pas avantage à choisir z 
-comme fonction inconnue de # et de O, la surface de contact étant 
caractérisée par la valeur constante @};, et les extrémités, par les fonc- 
tions données ,(f), @,(f) ; constantes dans l'exemple traité par 
NEUMANN ou par STEFAN. 
Ce changement de la variable indépendante et de la fonction 
-donne, sans difficulté, les équations non linéaires 


k, d'2 __ oz (EY 1 
C9 087 ot \ 28 | mis 
R 02 __ 02 nr . 
C…0 287 of | je 
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et à la surface de contact  — 0, 


Z(ê) — 2:(?) 
I 02 I 021 Lo 92, 22 04 
Ra 28: ki 20: Fa kik 29, 28 ot à 


La transformation qui fixe les limites des domaines de z, et de z,, ne 
facilite en rien la solution. 

Et si on présente à un mathématicien, ces nouvelles équations, il 
songe à les ramener à une forme linéraire, plus facile à étudier quand 
on ne se préoccupe pas des conditions limites, — et fournit les équa- 
tions mêmes qui m'ont servi de point de départ (I, $ 5)! 


14. — Mur homogène. — On appelle mur l’espace compris entre deux 
plans parallèles. Soit H la distance des deux plans. Supposons donnée 
la distribution initiale de température © (z, o) en fonction de z, dans 
l'intervalle O — H. Pour être maître des conditions frontières, je 
choisis la solution qui maintient la température zéro pour z — 0 et, 
aussi pour z == H : 


en posant 
Bo GYM) cp _(e+r+7H)% 
ty D= > Ce usa # ) 
; Nn = — © 
dont les termes se compensent deux à deux pour 2 — 0 et pour z — H. 
A cette intégrale, il faut ajouter celles qui proviennent des 
frontières 


n{ 2° 


e(z, t) = | 8(0, r)e : ak(t— 7) x ay cp dr 


F2 21rkË — T)"/2 
ni &—H}æ = 
* _ ak(t—5) (Zz — H)Vce dr 
ï A MU: DE QU 


qui donnent o pour l’état initial. 

I1 paraît inutile d'écrire les équations de problème de fusion ; 14° 
forme du noyau %, ne permettra évidemment pas de deviner, comme 
NEUMANN et STEFAN, la loi de propagation Ë(f) de la surface de fusion, 
même dans les cas les plus simples où @ (7. 0), @ (0, r), © (H, tr) seraient 
des constantes. 
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15. — Problème shhcrique tcnrope. — Amvonsan problème spèé- 
rique. qui, sul, présente avec le problème terretr assez d'analoes 
géométriques, pour guider, — si on savait le récemre — les vues 
qualitatives, et un peu quantitatives, ser ka marche de ln coma 
tion profonde de notre globe. 

Malheureusement tous les essais sont restés infrectee. 

Le problème qui smble le plus simple, a 
celui de NEUMANX. Dans un mikeu Kquide infni à tempétratere uni 
forme, en refroiit imtantantment me foi à une température 
inférieure à la température de congélation, et on maintient cette te 
pératare constante. Il se formera un petit noyean sode snlhèrique : on 
demande comment il se développe : comment croit avec le temps, le 
rayon de sa surface de congélation. 

S$ on essaie d'imiter ls solution @ NEUMANX, on se trouve asset 
arrèté, parce que là solution sphérique analogue est 


4 x AS ra 
as SES 
; ke dénominateur r gène toute ka suite. 


16 — Le problème sphérique kotrope utile, analogue am problème 
plan, est celui dans lequel là surface r — R d'une sphère d'a, avant 
à l'origine du temps une dxtrnibution de température otrope fsverable, 
serait portée brusquement à une température @, inèreure à là tems- 
pérature de sokification, là température ancentre étant maistenue 
constante @. 

On aura dans ke giace 
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avec 
x 0, _ }, 26: a 
(+ ki 27 Fa Se Lo: 
' : ONE _… Al DE or 
Il y aurait fusion, si e était positif. Pour que l'épaisseur solidifiée 


croisse avec le temps, il faut que les températures extrêmes Oo, Or et 
l’état initial soient choisis de telle sorte que = soit négatif. 


Dans le problème de fusion, les fonctions données sont 


OR, t) glace 
8,(Y, 0) eau 

ie fonctions cherchées : 
À (y, 0) glace 
8(R, f) eau 


R(t) figure dans les équations XX et XXI, d’une manière qui ôte 
tout espoir d'en obtenir l'expression par une voie simple. 


17. — La formation d'équations intégrales équivalentes — ana- 
logues à celles du problème général du mur plan $ 14 — est intéres- 
sante en elle-même. 

La masse est supposée initialement liquide. On refroidit la surface 
extérieure R, au-dessous du point de fusion ; et on cherche la marche 
de la surface de séparation du solide et du liquide versle centre. 

Les noyaux utiles sont un peu plus compliqués que pour le problème 


plan. 
Je pose 
k HR Q=r=nR aps — _(+r—nR) ap: 
JTE CAT LT 4kit et 4kit 
3 nn = — © 


Pour 7 — 0, N est évidemment nul, et linéaire en 7 pour 7 très petit. 
N Sete : es 
, st donc fini à l’origine, et convient au cas où il n’y a pas de source 


de chaleur au centre. 
A la surface libre, pour 7 = R, N est également nul, un terme # du 
premier groupe étant compensé par le terme (# + 2) du second groupe. 
N est le noyau qui correspond à l’état initial, donné seulement de 


OA R,. 
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Pour la température superficielle donnée, je pose 


k , 1 I ras T'aue* | 
XXI 1 Q/R Re 4hkité—T) _ akilt— 7 
( I) ee 1 e 
avec 2 termes seulement. 

M, est toujours nul au centre, et linéaire en r, au voisinage, ce qui 
laissera la température finie au centre. 


Mêmes formes Nr, Mn ( avec = pour l’eau centrale. 
18. — Je prends pour la glace (I) : 


rR nf 
Ge 4,97 [vit oNdr + 7 | RBUR, SM r 
(e] O 


La température superficielle est donnée ®, (R;f) ; mais comme la 
glace n'existait pas encore au temps £ — o, la température initiale 
Yi (r, 0) est une inconnue du problème. 

Pour l’eau (II) il y a permutation des fonctions connue et inconnue. 
La distribution initiale ® (7, o) est connue ; mais la température à la 
surface R, qui cesse d’être liquide dès le début, x (R, é) est une in- 
connue du problème. 

J'écris donc 


'R Vale) 4 
(XXV) er, =; je vb; O)Nndy + > J RU(R, *)Mnés. 
o o 


Les fonctions W;, Wir et R(f) devront être déterminées par ces deux 
équations XXIV, XXV, dont les 4 noyaux sont connus, et par les 
3 équations XX, XXI, à la surface de solidification 7 — & (é). 

La solution du problème plan n’a été devinée que dans le cas où 
l'épaisseur liquide initiale était infinie, ce qui réduisait à un seul terme 
chacun des noyaux équivalents à N et M. Mais si l'épaisseur d’eau à 
congeler était finie, la condition thermique à la 22 surface introduirait 
dans le problème plan un développement en série comme noyau, — et 
l'impuissance du mathématicien resterait la même que dans le pro- 
blème sphérique. 

Même dans le cas le plus simple : température initiale de l’eau On 
uniforme et température superficielle de la glace ®, constante, — la 
complication des noyaux, et l'ignorance sur les fonctions Wir, Y, 
nous laissent complètement désarmés. 
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19. — Peut être serait-il moins difficile d'attaquer les 5 équations 
d’une manière une peu différente, parce que les fonctions inconnues 
seraient moins engagées sous des signes algébriques. Ce serait de se 
donner arbitrairemet une loi de propagation de la solidification & (4), 
et de chercher à quel groupe d’états initiaux et de températures à la 
surface R, cette loi R (/) serait adaptée. On pourrait évidemment choi- 
sir, en outre, l’une des 4 fonctions ®, W ; par exemple prendre ®; (tem- 
pérature superficielle de la glace) constante. 


20. — La présence des deux fonctions inconnues W;, W;, semble une 
complication. 

On pourrait évidemment tirer des intégrales de la chaleur, une forme 
de solution satisfaisant directement aux équations XX, c’est-à-dire, 
fournissant l’état du liquide intérieur à une surface isotherme (@};) 
mobile (R (£)), et l’état de la glace entre la même surface isotherme mo- 
bile et la sphère R maintenue à température connue. Mais je n’ai réussi 
à formuler de cette manière rien qui mérite d’être transcrit. 


21. — Conclusion. — Voici donc un problème celui de la propaga- 
tion isotrope de la solidification par conductibilité, inattaquable par 
les méthodes que les mathématiciens ont déjà fournies aux physiciens. 

Toutes les difficultés proprement physiques du problème ont été 
supprimées. Il ne reste que le squelette de la question. Mais, le fait que 
le corps est limité, soit entre 2 plans, soit entre le centre et le rayon R, 
impose aux noyaux N, M, une forme compliquée qui interdit toute 
divination. 

Surtout, l’inconnue est le déplacement d’une frontière mobile. 
Aucune des méthodes en usage parmi les mathématiciens, ne s'adapte 
(même sans exigences de rigueur) à l’attaque de semblables problèmes. 

M. B. 


Note 1. — Pour toutes les formules utilisées dans cet article, — dues 
à FOURIER et LAPLACE, —'ou POISSON — il est commode de se reporter 
à un memento très bien établi par sir W. THOMSON vers 1850 pour son 
usage personnel : « Compendium of the Fourier mathematics for the 
conduction of heat in solids.. » Mat. Phys. Papers, t. II, p. 41 à 60; 
voir aussi p. 63-65. 
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On consultera avec intérêt : 

LAPLACE. — Mém. sur divers points d'analyse, Journal de l'Ec. Pol., 
XVE cahier (en particulier p. 235-244, et en genéral : Œuvres. 

FourtER. — Théorie de la Chaleur et Œuvres. 

Poisson. — Théorie Mathém. de la Chaleur et Œuvres. 

RIEMANN-WEBER (1901). — Die partiellen Differentialgleichungen 
der Mat. Phys., t. II. 

J'ai mis à contribution ces diverses sources pour constituer le Me- 
mento qui suit. 


A. — DÉBIT DE CHALEUR DONNÉ A LA FRONTIÈRE 


Plan. — Une source plane snstantanée, qui fournit à l’époque é{ = 0 
la quantité de chaleur g par unité de surface, à un solide primitive- 
ment à la température uniforme 9, — et cesse aussitôt — produit 
dans ce solide une distribution de température 


x°co 


e 4 Jo 
(0) — 0, = 9 —— =. 
(x, ?) 0 q 2Vr Bt 
Point. — Une source ponctuelle (à l’origine) instantanée, qui fournit 


la quantité totale de chaleur Q et cesse aussitôt, produit dans tout leso- 
lide primitivement à @, une distribution isotrope de températures 


co “ Br°/o (E)” À 


\ 


O(7,1) — 8 —=<€X 


Sphère. Une source sphérique (rayon R) instantanée, qui fournit 
la quantité totale de chaleur Q, et cesse aussitôt, produit dans tout le 
solide, primitivement à @,, une distribution isotrope 


(r — R}*ce __r+R}% 
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Grâce aux deux exponentielles, la température à l’origine (7 — o) 
reste finie et égale à 


nue PEER 
et SE 


SUR QUELQUES PROBLÈMES DE LA PHYSIQUE MATHÉMATIQUE 


B. — TEMPÉRATURE DONNÉE À LA FRONTIÈRE 


Plan, — Te plan x — 0, primitivement à la température o est porté 
brusquement à la température ©. 
La température était 0 partout. La température devient 


Renversons ia température O, au bout du temps 7 très court. Après 
ce renversement, la température à l’origine (x — 0) est redevenue zéro, et 
elle est partout 


x?cp 


D oct, ae se 


(forme obtenue d’une tout autre manière par W. THOMSON en 1855, 
TE 605: 

On tire facilement de là, que si, dans le milieu indéfini à o°, le plan 
x = X est porté aux températures O, (X, r) à partir du temps r = 0, 
la température d’un plan x quelconque, à l’époque £ est 


EX), ab) di 
O(x A) — NEtE ; e,(X te E—) 
C. — ETAT INITIAL DONNÉ 
Plans. — La solution est due à LAPLACE (Mém. sur divers points 
d'Analyse). 


7 7 
A _ bad 
: re Ze f ME à #4 Va)e ds 


L'état initial @ (x, 0) est F (x), comme on le vérifie facilement. On 


s'assure facilement aussi — en dérivant sous le signe | et intégrant par 


parties la dérivée par rapport au temps, — que l’équation aux dérivées 
partielles est satisfaite. 

Espace 1ndéfini du côté x > 0.— La température initiale n’est connue 
que du côté positif ; d’autre part, il est commode, pour le libre choix 


MARCET, BRITLLOUIN 


de la condition frontière, de rendre constamment nulle la température 


à l’origine. C’est ce qu’on obtient en prenant F (— x) = — F (x), 
| ,* /cp | 00 
D ni : V Le 2q VE)e “a 20 LÉ (+ se 4VÀ) Je * à 
Vr _œ C? x VE ( 


ou, ce qui revient au même : 


200 ü ‘ 2 
DES cp / | INT A CRT ù 
0 — NE 11 F(y, o)le e dy. 


Remarque. — Cette transformation de l’exposant, avec changement 
de la variable d'intégration (ou v. v.) est constamment utilisée pour 
reconnaître les valeurs limites des intégrales, dans la théorie de la 
chaleur. 

On passe sans difficulté de cette dernière forme au noyau N du mur 
homogène $ 14. 

Sphère. — La distribution isotrope ne fait connaître les températures. 
initiales F (‘, 9) que pour 7 => 0. On prendra donc, si l’espace a une éten- 
due infinie 


OP RS à | RENE See Sa 
8 — — NE É YE(y, 0) Le Le 4h dy 

La température reste finie à l’origine, où il n’a y pas de source de 
chaleur. 

Si la sphère est limitée par une surface de rayon R, l’espace intérieur 
(o < 7 < R) est comparable au mur plan. D'où la forme de série pour le 
noyau relatif à l’état initial ($ 17, $ 18). 


(Deux conférences données à l’Institut H. Poincaré le 13 et 
16 Mai 1929 et rédigées par M. Coulomb). 
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